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Abstract 

In the present paper we prove a form of Arnold diffusion. The main result says 
that for a " generic" perturbation of a nearly integrable system of arbitrary degrees of 
freedom n 2 

H (p) + eHi(9,p,t), 6eT\ p e B n , teT = R/T, 

with strictly convex Hq there exists an orbit (0 e ,p e )(t) exhibiting Arnold diffusion in 
the sens that 

Bu P ||p(i)-p(Q)||>j(zri)>o 

t>0 

where l{Hi) is a positive constant independant of e. 

Our proof is a combination of geometric and variational methods. We first build 
3-dimensional normally hyperbolic invariant cylinders of limited regularity, but of 
large size, extrapolating on |Be3] and [KZZj . Once these cylinders are constructed 
we use versions of Mather variational method developed in Bernard |Bel] , Cheng- Yan 
[CYHIU72] . 

1 Introduction 

Let (0,p) € T n x U be the phase space of an integrable Hamiltonian system Hq(p) with 
T n being the n-dimensional torus T n = W 1 /II 1 3 9 = (Ox, - ■ ■ ,0 n ) and U being an open 
set in M n , p = (pi, • • • ,p n ) £ B n . Assume that Hq is strictly convex, i.e. Hessian dp. p .Ho 
is strictly positive definite. 
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Consider a smooth time periodic perturbation 



H e (e,p,t) = H ( P ) + eH 1 (e, P ,t), iGT = R/T. 

We study Arnold diffusion for this system, namely, existence of orbits {(9,p)(t)}t such 
that 

\p(t) — p(0)\ > O(l) independently of e. 

Integer relations k\ ■ 8 p Hq + k$ = with k = (ki,ko) G (IT 1 \ 0) x Z and • being the 
inner product define one- dimensional resonances. Under the condition that Hessian of Hq 
is non-degenerate, these resonances define smooth hyper-surfaces embedded into action 
space 

r r {peB n : h ■ d p H + ko = 0}. 

If one intersects k linearly independent resonances {/cj}" =1 , we get a k-dimensional reso- 
nance nr^ , which is defined by an (n — /c)-dimensional surface in B n . 

1.1 Apriori unstable systems 

In the case n = 2 Arnold proposed the following example 

I 2 p 2 

H(I tp,q,p,t) = — + — + e(l -cost/)(l + /u(sin6> + sint)). 

The feature of this example is that it has a 3-dimensional normally hyperbolic invari- 
ant cylinder (NHIC). This allows in a sense to reduce 5-dimensional dynamics to a 3- 
dimensional one. There is a rich literature on Arnold example and we do not intend to 
give extensive list of references; we mention [AKN] IBB} IBe4l IBsll !Zha , and references 
therein. 

This example gave rise to a family of examples of systems of n + 1/2 degrees of freedom 
of the form 

H s (I,<P,P,q,t) = H (I) + K (p,q) + eH^I, p,p, q, t), 

where (q,p) £ TP -1 x M n_1 , I £ R, 93, t £ T. Moreover, the Hamiltonian Ko(p,q) has 
a saddle fixed point at the origin and Kq(0, q) attains its strict maximum at q = 0. For 
small e this system has a 3-dimensional NHIC A. 

For n = 2 systems of this type were successfully studied by different groups. Two 
groups were using deep geometric methods. 

- In [PHI IDGLSl IDLS] the authors carefully analyze two types of dynamics induced 
on the cylinder A. These two dynamics are given by so-called inner and outer maps. 

— In [Tl| IT2J a return (separatrix) map along invariant of A is constructed. A detailed 
analysis of this separatrix map gives diffusing orbits. 

The other two groups [Bell ICYlj are inspired and influenced by Mather variation 
method [Mai, Ma2] and build diffusing orbits variationally. We essentially rely on their 
technique in this paper. 
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In the case n > 2 Treshev [Tit IT2] and Cheng- Yan [CY2] also managed to prove 
diffusion. In Treschev's paper he even showed existence of orbits with an optimal diffusion 
speed ~ | lne|/e. Other examples of Arnold diffusion can be found in |Bs2} IBo( IBK[ iKLTj 

iin^iKmii^ii^raiM^imizh^] . 

1.2 Dynamics along a single resonance for n = 2 and crumpled invariant 
cylinders 

Pick a single resonant line C B 2 . Notice that on a dense set in there is an additional 
resonant relation. If one deletes a finite number of such additional resonant relations 
with relatively small fc's usually called double resonances, then on each segment between 
consecutive deleted points one can hope to find a "nice" smooth NHIC. This turns out to 
be wrong] However, one can indeed find a NHIC whose regularity blows up as e — > 0. 
Moreover, one can use this cylinder for Arnold diffusion. This is a new phenomenon 
discovered in this paper. 

1.3 Dynamics along a codimension one resonance for n > 2 and crumpled 
invariant cylinders 

Pick a codimension one (or dimension n— 1) resonant line T = n^Z^T^ C B n with {kj}jZl 
being linearly independent. As before on a dense set in V there is an additional resonant 
relation. However, qualitatively the picture as in the case n = 2. Namely, if one deletes a 
finite number of such additional resonant relations with relatively small fc's, then on each 
segment of T between consecutive deleted points there is a crumpled 3- dimensional NHIC 
A. It is a crumpled in the sense that its regularity blows up as e — > 0. With some efforts 
this allows to reduce dynamics to 3-dimensional one and essentially reduce the proof to 
the twist maps case. 

1.4 Main result 

We study dynamics near a resonance of codimension one, i.e. near a segment in B n . For 
any resonance of codimension one there is an integer linear symplectic transformation 
which brings integer vectors k±, . . . , k n _\ E Z n , defining the resonance, to the form kj = 
(0, • • • , lj, 0, • • • ,0). Since we are interested in a local property assume that a resonance, 
denoted T, of codimension one is of the following form: 

(d Pl H (p),--- ,d Pn _ 1 H (p)) = (e 1 ,--- ,9 n -i) = for e = 0. 

In the case Hq(j>) = | Y^=iPj we have T = {(pi, • • • ,p n -i) = 0}. Thus, it is naturally 
parametrized by p n . 

Consider the space of C r perturbations C r (T n x B n x T,]R) with a natural C r norm 
given maximum of all partial derivatives of order up to r. Denote by S r the unit sphere 
in this space. 
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Theorem 1.1. For r ^ 4 ; there is an open and dense setlA C S r , a nonnegative function 
I : S r — > M + with l\u > and a positive function e$ = Eq{Hi), we write V = {eH\ : H\ € 
U, < e < £$}. We have that, for an open and dense set of eH\ £ V the Hamiltonian 
system H t = H$ + eH\ has an orbit {(6,p)(t)} t whose action component 

|b(T)-K0)|| 

Moreover, for all < t < T the action component p(t) stays close to the codimension one 
resonance T. 

Remark 1.1. This Theorem provides a form of Arnold diffusion for generic Hamiltonian 
systems. The type of generic condition in Theorem ll.il is a version of Mather's cusp 
residue condition introduced in JMa3l. 

The present work is in large part inspired by the work of Mather |Ma3l IMa41 IMa5j . 
In [Ma3|, Mather announced a much stronger version of Arnold diffusion for n = 2 (the 
system is time-periodic hence the degree of freedom is 2^). The proof of Mather's result 
is partially written (see |Ma4] ) . and he has given lectures about some parts of the proof 
[Ma5j . One of the ideas underlying his proof is to construct diffusion along a segment of a 
resonance and away from other low order resonances. Conceptually the proof of our result 
has similar features to parts of Mather's proof [Ma4j for single resonances. The novelty 
of our approach is the use of normal form theory and construction of normally hyperbolic 
cylinders in an a priori stable setting. Application of normal forms to construct normally 
3-dimensional hyperbolic invariant cylinders in apriori stable situation in 3 degrees of 
freedom is proposed in [KZZJ. Independently in the case of arbitrary degrees of freedom it 
is proposed in |Be3j . In the latter it is shown that such cylinders have length independent 
of e. 

1.5 Plan of the proof 

The proof of this Theorem proceeds in three steps. 

Step 1. Build a normal form for H £ for for p near T. In section [3] we prove the existence 
of a normal form, which takes a particular nice form along subsegments of V, which we will 
call passage segments, defined in the next section. The length and choice of the passage 
segments depends on Hq and H\ only. 

Step 2. For H\ 6 W, we establish existence of finitely many 3-dimensional normally 
hyperbolic cylinder along V. This is discussed in Section 01 

Step 3. For a generic perturbation, we show that there exists diffusion orbit along 
a passage segment, using the normally hyperbolic cylinders. This steps uses variational 
methods of Bernard [Bel] and of Cheng- Yan |CY1[ ICY2j which are based on ideas of 
Mather (see |Ma4] ) . These constructions are discussed in Section [5] and Section [6j 
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2 Notations and terminology 

We denote 9 s = (6>i, • • • , 6 n -i), P s = (pi, • • • ,Pn-i) and 0-* = 6* n , pf = p n . These are the 
slow-fast variables associated to the resonance V = {d p sHo(p) = 0}. It is natural to use pf 
as a parameter for T, i.e. we may write Tf]B = {p*(p f ) = (pt(p f ) , P f ) , P f G [a m i n , a max ]} . 

Z(9 s ,p) := JJ H 1 (9 s ,p s ,9 f ,p f ,t)d9 f dt. 

If the perturbation H\{Q,p, t) is expanded as 

ffi(0, Pl t) = H x {e s ,e f ,p,t) = Yl h [ ks,kf,i ] (p)e 2ln{ka - es+k, - 9f+l - t \ 

jfc»ez r »- 1 ,jfc/ez,Je2 

then 

z(r,p) = ^/ l[fc , i0i0] ( P )e 2 -( fcfl -^. 

We would like to impose the following set of non-degeneracies and notations. Consider 
the function Z(9 S ,p*(pf)) as a family of functions on T™^ 1 parametrized by pf . 

Call a value pf on T regular if Z(9 S ,p*(pf)) has a unique global maximum on T s 9 9 s 
at some (9| = 9 s (pf ). We say the maximum is non-degenerate if the Hessian of Z with 
respect to 9 s is strictly negative definite. 

Call a value p> on T bifurcation if Z(9 s ,p*(pf)) has exactly two global maxima on 
T s 3 9 s at some 0f = 0f(/) and 6>| = 6»|(/). 

Call a regular on T non- degenerate if the unique maximum is non-degenerate. If p? is 
a bifurcation, it is called non-degenerate if both maxima are non-degenerate, furthermore, 
the values at these maxima moves with different speed with respect to the parameter pf 
Otherwise, it is called degenerate. 

The generic condition that defines U C S r is a higher dimensional version of the con- 
ditions (C1)-(C3) given by Mather [ Ma3j . These conditions may be described as follows: 
Each value G [a m in, Umax] is a non-degenerate regular or bifurcation point. Note that 
the non-degeneracy condition implies that there are at most finitely many bifurcation 
points. Let a± < ■ ■ ■ < a s -\ be the set of bifurcation points in the interval (a m i n ,a ma x), 
and consider the partition of the interval [a m i n ,a max ] by {[aj, aj+i]}*^. Here we give an 
explicit quantitative version of the above condition: There exists A > such that 

[GO] There are smooth functions 9j(p>) : [aj — A, aj+i + A] — > T n , j = 0, • • • , s — 1, such 
that for each pf G [aj — A, dj+i + A], 9j(pf) is a local maximum of Z(9 S ,p*(p*)) 
satisfying 

XI ^ -dUsZ{9 s j ,p) ^ J, 
where I is the identity matrix. 

[Gl] For pf G (aj,aj + i), 9? is the unique maximum for Z. For pf = aj + \, 0? and 0j +1 
are the only maxima. 
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[G2] At pf = Oj-fi the maximum value of Z has different derivatives with respect to pf , 
i.e. 

^Z(0*(a i+1 ),p*(/)) + -^Z(9 s j+l (a j+l ),p*(pf)). 

Theorem 2.1. The set IA of functions H\ € S r such that the corresponding Z(9 s ,p) 
satisfies conditions [G0]-[G2] is open and dense. 

The proof of Theorem 12.11 will be given in Appendix [Al 

Write oj{p) = d p Ho(p) = (d p sHo,d p fHo), clearly for any p £ T we have that u(p) = 
(Q,d p fHo)). We say that p$ has an additional resonance if there exists integers k n ,l such 
that k n d f Hq(p) + 1 = 0. Given a large integer K, let 

z K = { P ernB; 3k n ,i e z, \k n \, \l\ ^ K, k n ■ d pf H ( P ) + l = o}. (l) 

Given Hi £ U, we will define a small 5 = 5(Hi,n,r) > and integer K = K(5,n,r) and 
call the elements of T,k punctures. We need to exclude a neighborhood of the punctures 
from rnB. Let U i stand for 3ee neighborhood of then F (1 B \U i 

is a collection of disjoint segments. Each of these segments is called a passage segment. 
On a neighborhood of each passage segment there exists a convenient normal form for the 
Hamiltonian H e . 



3 Normal forms 

Let T = {(p s = p*(pf))} be the resonant line of equation 8 p sHq = 0. For pGTwe have 
u{p) = (0,d p fHo). We say that p has an additional resonance if the remaining frequency 
d p fHo(p) is rational. In order to reduce the system to an appropriate normal form, we 
must remove some additional resonances. More precisely, let T>(K, s) C B be the set of 
momenta p such that 

• Hdps-ffoWH ^ s, and 

• \k f d pf H (p) + k*\ ^ 3Ks for each (k f , k l ) G I? satisfying max(V , k l ) e]0, tf]. 

Theorem 3.1. [Normal Form] Let Hq(p) be a C 4 Hamiltonian. For each 5 G]0, 1[, there 
exists positive parameters Kq,€o,(3 such that, for each C 4 Hamiltonian Hi with \\Hi\\ C 4 ^ 
1 and each K ^ Ko,e ^ eo, there exists a C 2 change of coordinates 

$ : T" x B x T — ^fxfxT 

satisfying \\<& — id\\co ^ ^fk and \\ & — id\\(j2 ^ 5 and such that, in the new coordinates, the 
Hamiltonian Hq + eHi takes the form 

N e = H (p) + eZ(e s ,p) + eR(6,p,t), (2) 

with \\R\\ C 2 O on T n x V{K,fi<i l l A ) x T. We can take K = c5~ 2 ,f3 = c^'Vo = 
<5 6n+5 /c, where c > is some constant depending only on n and \\Hq\\ C 4,. 
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The proof actually builds a symplectic diffeomorphism 6 of T n+1 x R n+1 of the form 
$(8,p,t,e) = {<S>(9,p,t),e + f(9,p,t)) 

and such that 

N e + e = {H t + e)o <j>. 
We have the estimates ||$ — z^Hco ^ \fe and ||$ — irf||c2 ^ 5. 

Remark 3.1. [Length of passage segment] On the interval, the distance between 2 adjacent 
rationals with denominator at most K is It follows that the distance between p{,p{, E 

Y<k (see ([JJ J is at least \\d 2 H ' 1 \\^ ^ c(c(n,r)^ 1 S)' F:: ^ , assuming that \\d 2 HQ 1 \\ is bounded 
by some universal constant. 

To prove Theorem 13. II we proceed in 3 steps. We first obtain a global normal form N e 
adapted to all resonances. We then show that this normal form takes the desired form on 
the domain T>k,e- However, the averaging procedure lowers smoothness, in particular, the 
technique requires the smoothness r ^ n + 5. To obtain a result that does not require this 
relation between r and n, we use a smooth approximation trick that goes back to Moser. 

3.1 A global normal form adapted to all resonances. 

We first state a result for autonomous systems. The time periodic version will come as 
a corollary. Consider the Hamiltonian H e ((p,J) = Hq(J) + eHi(ip,J), where (</?, J) E 
T m x R m (later, we will take m = n + 1). Let B = {\ J\ ^ 1} be the unit ball in R m . Given 
any integer vector k E Z m \ {0}, let [k] = max{/cj}. To avoid zero denominators in some 
calculations, we make the unusual convention that [(0, • • • ,0)] = 1. We fix once and for 
all a bump function p : M. — > R be a C°° such that 

p(x) = 

and < p(x) < 1 in between. For each (3 > and k E Z m , we define the function 
Pk(J) = P(. ^j'S'/^[k] )' wnere /? > is a parameter. 

Theorem 3.2. There exists a constant c m > 0, which depends only on m, such that the 
following holds. Given : 

• A C 4 Hamiltonian Hq(J), 

• A C r Hamiltonian Hi(ip,J) with \\Hi\\c- = 1, 

• Parameters r ^ m + 4, 5 e]0, l[, e e]0, 1[, /3 > 0, K > 0, 
satisfying 
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-1 

• K ^ c m 5 r - m ~ 3 , 

i/iere exists a C 2 symplectic diffeomorphism such that, in the new coordinates, the 
Hamiltonian H t = Hq + eHi takes the form 

H e o<5> = H + eRx + eR 2 

with 

• R± = ^2 ke i m \k\<x Pk{J)hk(J) e2 k i here /ifc(J) is the k th coefficient for the Fourier 
expansion of H\ , 

• H-^llc ^ 8, 

• ||$ - id\\ c o < 5y/e and ||$ - id\\ C 2 ^ 8. 

We now prove Theorem 13.21 To avoid cumbersome notations, we will denote by c m 
various different constants depending only on the dimension m. We have the following 
basic estimates about the Fourier series of a function g((p,J), Given a multi-index a = 
(ax, • • • , a m ), we denote \a\ = ax + • • • + a m . 

Lemma 3.1. For g(<p,J) 6 C r (T m x B), we have 

1. Ifl^ r, we have \\gk(J)e 2lTi{k ^\\ c i < [Jfe]'- r ||s||c- 

2. Letgk(J) be a series of functions such that the inequality \\djagf.\\ c o ^ M[/c]~l a l _m_1 
holds for each mult-index a with \a\ ^ I, for some M > 0. Then, we have 

II J2 k ez^9k(J)e 2 ^ k ^\y ^ cK m M . 

3. Let U+g = Y,\k\>K 9k{J)e 2lTi{k ' v) ■ Then for I ^ r - m - I, we have ||n+g|| ci 

v T^m— r+l+X II „|| 
K m K ^ \\9\\C r - 

Proof. 1. Let us assume that fc ^ and take j such that kj = [k]. Let a and r\ be two 
multi-indices such that \a + r/\ ^ I. Finally, let b = r — I, and let f3 be the multi-index 
/3 = (0, . . . , 0, b, 0, . . . , 0), where = b. We have 

g k (J)e 2in(k '^ = [ g(6,J)e 2i7r ^- e) d9= [ g(6 + <p, J)e- 2i7r(k ' 8) dd, 

hence 

d^jv( 9 k(J)e 2 ™^) = I d^j,g(d + v,J)e- 2 ^ k ' e Ud, 



T" 



T» 



9^^9(9 + <p, J)_ e -2vK{k,e) de 



(2i7rkj ) 
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Since \a + /3 + r]\ ^ r, we conclude that 

\\g k (J)e 2i < k ^\\ cl < \\g\\cr/(Mk]) h < \\g\\ C r[k] l - r . 

2. We have \\g k { J) e 2ln( - k ^ \\ c i < 

|| h k {J)e 2 ^ k ^\\ c i ^ J2 ci\k\~ r+l M ^ Cl K m M, 
fcez m fcez m 



m—1 

uiiau n, m — Z^keZ m l a l 

3. 



recall that K m = X^fcez™ W 



l|n+5lb^c E \k\- r+2 \\g\\ C r^cK- r+m + 3 Yl W^Mcr 
\k\>K \k\>K 

^ cK~ r+m+3 Km \\g\\ C r = c Km K- r+m+3 \\g\\ C r. 

3. Using 1., we get 

iin^ib < E { k t r b\\cr < ii 5 iic^ m ' r+i+1 E 

|fc|>K |Jfc[>A- 

^ ii 5 iic^ m ~ r+ ' +1 E M -m_1 - 

□ 

Proof of Theorem VJJA Let G(v?, J) be the function that solves the cohomological equation 

{Ho, G} + Hi = R\ + R + , 

where -R+ = ILtili. Observing that pk{J) = 1 when k • 8jHq = 0, we have the following 
explicit formula for G: 



n(,„ t\ (1 ~ Pk(J))hk(J) 2m(k-tp) 



where each of the functions (1 — pk(J))hk(J)/{k ■ 8jHq) is extended by continuity at the 
points where the denominator vanishes. This function hence takes the value zero at these 
points. G is well defined thanks to the smoothing terms 1 — pk we introduced, as whenever 
k ■ 8jHq = we also have 1 — p k = and that term is considered non-present. 

Let <£* be the Hamiltonian flow generated by eG. Setting Ft = R\ + R+ + t(H\ — R\ — 
we have the standard computation 

d t ((H + eF t ) o $*)) = ed t F t o $* + e{H + eF t , G} o 

= e(d t F t + {# , G}) o $* + e 2 {F t , G} o $' 
= e 2 {F t ,G}o<D*, 
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from which follows that 



H e o = H + eRi + eR + + e 2 f {F t , G} o $*di. 

J o 

Let us estimate the C 2 norm of the function R2 := -R+ + e J 1 {i 7 t,G'} o <E>*<ii. It follows 
from Lemma 13. II that 



\R+\\ C 2 < K mJ FC- r+m+2 ||F 1 || c , < -6. 



We now focus on the term L{Ft, G} o To estimate the norm of Ft, it is convenient 

to write Ft = Ft + (1 — t)Ri, where Ft = (1 — t)R+ + tffi. Notice that the coefficients of 
the Fourier expansion of Ft is simply a constant times that of H±, Lemma [3. II then implies 
that 

Halloa < W 3_r H^lll^ = «m||-Hl||cr- 
kez m 

provided that r ^ m + 4, where K m = Yl,z n W\ rn+l - 

We now have to estimate the norm of R\ and G. These estimates require additional 
estimates of the smoothing terms pk as well as the small denominators k-djHo. We always 
assume that I £ {0, 1, 2, 3} in the following estimates: 

- Pk (J) ± 1 =► I (A; • djH y l \ sC ^e-^lfcl" 1 . 

- [|(fc • SjiZoMc < c m j3- l - l e- ( - l+1 ^\\H \\ l + 1 on {p, ^ 1}. 

- \\p k (J)\\& < c m r'e^ /4 ||^o|lc4 and ||1 - p fc (J)|| C j ^ c^e-^H^oH^- 

We have been using the following estimates on the derivative of composition of functions: 
For f : — y k and 5 : R m — > M. m we have ||/ o 5 || c , < c m j\\f\\ c i (1 + Nl^)- 

- For each multi-index \a\ ^ 3, we have that 

\\d Ja ((1 - Pk(J))h k (J)(k ■ djHor 1 ) ||c° 

< I' 1 ~~ Pk(J)\\c\"i\ ll^llcl«2l \\( k • ^- Hr o)~ 1 ||cl»3l({ Pfc ^l}) 

ai+ot2+a:j=a 

ai+at2+a3=a 

. ) g-l«3l-i c -(l«3l+i)/4||H ||g»l+ 1 ) 
<c m/ 9-l a l- 1 €-f° [ l +1 V 4 [A ; ]l a l- r [|^o||^! +1 ||^i||cr. 

In these computations, we have used the hypothesis /3e 1//4 ^ ||i/o||c 4 - Since G((p,J) = 
T.kez^i 1 ~ Pk(J))hk(J)(k ■ djHo^e 2 ™^, Lemma ED implies (since r ^ m + 1) : 
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- \\G\\ cl s: c^r'-^-^^ll^oll^ll^illcr < e" 1 - 

We now turn our attention to R± = X^/tlcK pk{J)hk{J)e 2%1T ^ k ' ip ' > : 

- \\h k \\c < [^'"ll^illc- 

- \\Pkh k \\ cl ^c m ^e-'/ 4 [A;]- r+ '||iI Q |^ 4 ||iTi||c- 

- ||-Ri||c/ < c m f3~ l e~ l / 4: \\Ho\\ l ci \\Hi\\cr , provided r ^ m + 4. 
We obtain 

ll-^llc' ^ ll-^i||c ! + ll-^llc ^ c mfi~~ l e~ l ^\\HQ\\ l ci \\Hi\\c-r, 

and 

IK-Fis C}|lc 2 ^ ^ ll-Pi llci Q ii l|C|lci Q 2l ^ c m/5 4 e 1 ||i?o||c ,4 ll-^illc'''- 

|ai+c«2|^3 

Concerning the flow we observe that ||eG||^3 ^ 1, and get the following estimate (see 
e. g. PHI): 

- \\&-id\\ C 2 < c m e||G|| C 3 < c m /T 4 1 1 Fo|| 4 4 | \H X \\ C r ^ 6, 

- ||$* - id||oo < c m e||G|| c i < c m /3~ 2 Ve||Fo|lc4ll^illc 2 < <ty*- 
Finally, we obtain 

e\\{F t ,G}o&\\ c2 ^ Cm e\\{F t ,G}\\ c2 \\$% 2 

^ C m /3~ 4 1 1 Hq 1 1 ^4 || Ifl || £r < 5/2. 

□ 

3.2 Normal form away from additional resonances 

We now return to our non-autonomous system and apply Theorem 13.21 around the reso- 
nance under study. To the non-autonomous Hamiltonian 

H e (6,p,t) = H (p) + eH 1 (9,p,t) : T" x R" x T — > R 

we associate the autonomous Hamiltonian 

He(<P,J) =H (I) + e + eH 1 (9,I,t) : T n+1 xf +1 — ► R, 

where ip = (6, t) and J = (I, e). We denote the frequencies ui G M n+1 by ui = (a)*, oj s , ui 1 ) G 
W 1 " 1 xixR, and define the set 

Q(K, s) := {to G E n+1 : \\uj s \\ > s, \k f u f + fcV| > 3sK V(fc s , fc*) G Z^}, 

where we have denoted by the set of pairs (kf , k l ) of integers such that < max(&;^, k t ) ^ 
K. Note that 

V(K, s) = {p€R n : (d p H (p), 1) G a)}. 
11 



Corollary 3.2. There exists a constant c n > 0, which depends only on n, such that the 
following holds. Given : 

• A C 4 Hamiltonian Hq(p), 

• A C r Hamiltonian Hi(8,p,t) with = 1> 

• Parameters r ^ n + 5, 5 G]0, 1[, e G]0, 1[, /3 > 0, if > 0, 
satisfying 

-i 

• -KT c n 5 r -™- 4 , 



/3e 1/4 < ||#o I 



there exists a C 2 symplectic diffeomorphism <I> of T n+1 x M n+1 suc/j i/iai, in the new 
coordinates, the Hamiltonian H e = Hq + eH\ takes the form 

N e = H + eZ + eR 2 , 

with 

• H-Rallc- <S onT n x V(K,f3e 1 / 4 ) x T, 

• ||$ - id\\ c o ^ 5y/e and ||$ - irf||c<2 ^ <5. 
JTie diffeomorphism $ is o/ i/ie forme 

m P ,t,e) = ($(0,p,t),e + f(6,p,t)) 

where is a diffeomorphism ofT n x M. n x T fixing the last variable t. 

Proof. We apply Theorem 13.21 with H t , m = n + 1 and 5 = 5/2. We get a diffeomorphism 
$ of T" +1 x M n+1 as time-one flow of the Hamiltonian G. By inspection in the proof 
of Theorem 13,21 we observe that G does not depend on e, which implies that $ has the 
desired form. We have 

H £ o $ = flo(J) + efli + e«2 

where ||i?2||c 2 ^ <V^ and 



Ri(9,P,t)= P\ 



kf • dpfHp + k" dp. H + >fc . (Sit) 



[fe] 5 

Let us compute this sum under the assumption that p G *D(K, /Se 1 / 4 ). We have 

fc^ • d pf H 
I3e l ^[k] 



^ 1 
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hence 



V -dpfHo + k'dpsHo + k 1 



for k such that k s = = k t . For the other terms, we have, by definition of Q(K, s) 



k s d p sHo + k l 



k s d p sHQ + k l 



hence 



kf ■d pf H + k a d p sH + k t 



> 2 



/3e 1 /4[fc] 

and these terms vanish in the expansion of K\ . We conclude that 



Ri{9,p,t)= 9(k f ,o,o){p) 



hence R\ = Z — 11^ (Z), with the notation of Lemma 13. II Finally H e o <E> = Hq + eZ + ei?2 
with R2 = R2 — ITft-Z. From Lemma l3.lt we see that 



\\U+Z\\c2 ^ c n K m+3 - r \\Z\\ C r ^ CnK^^WH.Wcr ^ c n K m+3 - r ^ 5/2. 
On the other hand, ||i?2||c 2 ^ "V^j hence ||i?2||c 2 ^ 
3.3 Smooth approximation 

Finally we remove the restriction on r by the following smooth approximation lemma: 



□ 



Lemma 3.3. JM/ Let f : R n — > R be a C r function, with r ^ 4. Then for each r > 
there exists an analytic function S T f such that 

\\S T f-f\\ c3 <c(n,r)\\f\\ c3 r r ~ 3 , 

\\S T f\\ C ri<c(n,r)\\f\\ cl T-^- r \ 
for each r\ > r, where c(n, r) is a constant which depends only on n and r. 

If r < n + 5, we use Lemma 13.31 to approximate H\ by an analytic function H^. We 
can then apply Corollary 13.21 to the Hamiltonian 

H* := Hq + eH\ = Hq + £2^2 

with H2 = Hl/\\H\ \\c r 2 , with €2 = ^\\H\ \\c r 2 , and with some parameters r2 ^ r and 
82 ^ 5 to be specified later. We find a change of coordinates $ such that 

H* o $ = H + e 2 Z 2 + £2-^2 
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and ||i?2||c 2 ^ ^2 j where Z 2 (9 s ,p) = J H 2 d0^dt. As usual, we have denoted by H* and 
Hq the automomized Hamiltonians H* = H* + e and Hq = Hq + e. With the same map 
we obtain 

H e o $ = H + eZ + eR 

with 

/? = \\HH\cr2R2 + (Z- Z*) + (Hi - H^ o 

In the expression above, the map <I> is the trace on the (9,p,t) variables of the map <£. 
Choosing r = 8\ 3 \ we get 



Ht-H^cs ^c(n,r 2 )5^ 3 
H{\\ C r 2 ^c(n,r 2 )5^ 3 

r-3 

Z* - Z\\ C 2 < \\Ht - #i|| C 2 s$ c(n,r 2 )5; 2 - 3 
*||oa < <*2 ^ <5 < 1, 



and finally 



We now set 



|-R||c 2 ^ c(n,r 2 )5. 



r-3 

r 2 -3 



62 = 5 r-a /c(n,r 2 ) ^ 5 



and get ||-R||c 2 ^ 8. To apply Corollary [32] as we just did, we need the following conditions 
to hold on the parameters: 



- K ^ c(n,r 2 )5^- 3 ^-"-^ , which implies K ^ c n 5^ n -\ 

- (3 > c(n,r 2 )(l + \\H q \\ c a)5~^) which implies f3 > c n (l + \\H \\ c a)S~ 1/2 , 

- fie 1 / 4 WHoWc^S 1 ^^ which implies fie 2 ^ < ||#o||c 4 - 

We apply the above discussion with r 2 = 2n + 5 and get Theorem 13.11 Note the estimate 



r 2- 



□ 
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4 Normally hyperbolic cylinders 



In this section, we study the Hamiltonian in normal form 

N e (9,p,t) = H (p) + eZ(6 s ,p) + eR(6,p,t). 

We denote as above by p%(p*) E M n_1 the solution of the equation d p s Ho(pl(p j) , p j) = 0. 
We recall also the notation p*(pf) := (p'l(Pf),pf ) from Section [2J Fixing parameters 

A e]0, 1], a~ < a + , 

we assume that there exists, for each pf E [a_ — A, a + + A], a local minimum 9 s (pf ) of the 
map 9 s i — > Z(9 S ,p*(pf)), and that 91 is a C 2 function of pf . We assume in addition that 

\I^d 2 eses Z(6t(p f ),p*(pf))^I 

for each pf E [o_ — A, a + + A], where as before / is the identity matrix. We shall at some 
occasions lift the map 91 to a C 2 map taking values in R n_1 without changing its name. We 
assume that ||Z!c 3 ^ 1; an d set ||-R||c 2 = Finally, we assume that D~ X I ^ d 2 p Ho ^ D I 
for some D ^ 1. To simplify notations, we will be using the O(-) notation, where / = 0(g) 
means |/| ^ Cg for a constant C independent of e, A, 5, n and r. In particular, we will 
not be keeping track of the parameter D, which is considered fixed throughout the paper. 

Theorem 4.1. There exists eo E]0, 1[ such that, if 

< e < e A 7/2 , < 6 < \^ X 2 , 

then there exists a C 1 map 

(e s ,P s )(9 f ,p f ,t) : T x [a- - A/2,a+ + A/2] x T — >■ T n_1 x M n_1 

such that the cylinder 

C = {(9 s , p s ) = (& s p P?)(9f,pf,t)); pf E [o_ - A/2,a+ + A/2],(^,i) E T x T} 

is weakly invariant with respect to N t in the sense that the Hamiltonian vector field is 
tangent to C. The cylinder C is contained in the set 

V :={(9,p,t);pf E [o_ - A/2, a+ + A/2], 

\\(9 s -9t(pf)\\^0(el /4 X), ||/-pJ(/)|KO(4 /4 A 5 /V/ 2 )}, 
and it contains all the full orbits of N e contained in V . We have the estimates 
\\Q s (9f,pf,t) - 9i(pf)\\ < 0(X- 1 8 + A- 3 / 4 ^), 
\\P s (ef,pf,t)-pl(pf)\\ < ^(A^ + A-^Ve), 



o 



WJ) 



0(X- 2 JI+X-^V5). 
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The proof of Theorem 14.11 occupies the rest of the section. The Hamiltonian flow 
admits the following equation of motion : 

6 = dps Hq + edps Z + edps R 
p s = —edgsZ — edesR 
< 0f = d pf H + edpf Z + edpf R . (3) 
pf = -ed 9 fR 
i = 1 

V 

It is convenient in the sequel to lift the angular variables to real variables and to consider 
the above system as defined on x ]R n_1 x 1 x R x I. We will see this system as a 

perturbation of the model system 

9 s = dpsH , p s = -ed 9 sZ , 6 f = d pf H ,pf = , i = l. (4) 

The graph of the map 

(0f, p f,t)^(et( Pf ),pi(p f )) 

on K x J x M is obviously invariant for the model flow. For each fixed pf, the point 
($*(P/)>P*(P/)) i s a hyperbolic fixed point of the partial system 

6 s = d p sH (p s ,pf) , f = -ed 6 sZ{0 s ,p s J) 

where p* is seen as a parameter. This hyperbolicity is the key property we will use, 
through the theory of normally hyperbolic invariant manifolds. It is not obvious to apply 
this theory here because the model system itself depends on e, and because we have to 
deal with the problem of non-invariant boundaries. We will however manage to apply the 
quantitative version exposed in Appendix iBl 

We perform some changes of coordinates in order to put the system in the framework 
of Appendix [Bj These coordinates appear naturally from the study of the model system 
as follows. We set 

B{pf) := d 2 ps p sH {p*{pf)) , A(p f ) := -$. fl .Z(0;(pO.P.(P / ))- 

If we fix the variable p* and consider the model system in (6 s , p s ), we observed that this 
system has a hyperbolic fixed point at (O s (p^),pl(p^)). The linearized system at this point 
is 

9 s = B{pf)p s , p s = eA(p f )9 s . 
To put this system under a simpler form, it is useful to consider the matrix 

L (pf) ■.= (b^VJCb^VJ^b^V))" 172 ^ 2 ^)) 172 
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which is symetric, positive definite, and satisfies L 2 (pf)A(pf)L 2 (pf) = B(pf), as can be 
checked by a direct computation. We finally introduce the symmetric positive definite 
matrix 

A(pf ) ■- L(p f )A(p f )L{p f ) = L' 1 {p f )B(p f )L- 1 (p f ). 
In the new variables 

t = L~ 1 {pf)e s + t- 1 ' 2 L{pf)p s , n = L~\pf)6 s -e- 1 / 2 L(pf)p s , 
the linearized system is reduced to the following block-diagonal form: 

i = e l / 2 k{pf)i , r, = ~e 1/2 A(p f )v, 

sec [Bc3 for more details. This motivates us to introduce the following set of new coordi- 
nates for our full system: 

x = L-\pf){6 s - ei(pf)) + e-WLtfw - pl(pf)) 
y = L~V)(0 S " 0*V)) " z- l,2 L(pf)(p s -pl(pf)), 
I = e-Wpf , 6 = jO f , 
where 7 is a parameter which will be taken later equal to b 1 ! 2 . Note that 

,1/2 



0" = 9t(e 1 / 2 I) + h(e 1 / 2 I)(x + y), 



V s =V^I 2 I) + *—L-\eV 2 I){ X -y). 



Lemma 4.1. We have A(p^) ^ \J XjD I for each pf G [a_,a + ]. 

Proof. The matrix A is symmetric, hence it satisfies A A*/, where A* > is its smallest 
eigenvalue. The real number A* is then an eigenvalue of the matrix 



similar to 



B 
A 



which is 

-A 

. Since both A and B are square matrices of equal size, we conclude 



that A„, 2 is an eigenvalue of A 1 B 1 . Since \\A 1 || ^ A 1 and \\B 1 || ^ D, we have 
A" 2 ^ IIA- 1 ^- 1 !! ^ DX^ 1 . We conclude that A* ^ y/X/D. □ 
The links between the various parameters e, 5, 7, A, p which appear in the computations 
below will be specified later. We will however assume from the beginning that 

S < ps^X , Ve^p , < 7 < A- 

Let us first collect some estimates that will be useful to see that the system ([3]) is 
indeed a perturbation of the model system. 



Lemma 4.2. On the domain \\x\\ ^ p,\\y\\ ^ p, we have the estimates 

-In _ im f£ || ^ 0(A~ 3/2 ), ||9 p /L- 3 / 2 || < 0(A- 3 / 4 ) 

r-ei\\ ^ oix-^p), \\p s -p s j < o{e l ' 2 P ). 



\\L\\ = 0(X~ 1/A ), \\L~ l \\=0(l), 

I^IK^' 1 ). \\pl\\c = o(i) 



17 



Proof. We recall that L = ( y B l / 2 (B l / 2 AB 1 l 2 )~ 1 l 2 B 1 / 2 )'. Since D~ l I ^ B < DI and 
A/ ^ A ^ /, we obtain that ||L|| ^ 0(A -1 / 4 ) and that ^ 0(1); using the expres- 

sion L~ l = (B- x l 2 (B x l 2 AB x l 2 ) x l 2 B- x l 2 ) Xj2 , we obtain that ||L|| < 0(A -1 / 4 ) and that 
^ 0(1), To estimate the derivative of L, we consider the map F : M i — > M 1 / 2 
defined on positive symmetric matrices. It is known that 

dF M -N = / e~ tMl/2 Ne- tMl/2 dt. 
Jo 

To verify this one can diagonalize M, perform integration, and match terms in (M 1 / 2 + 
edF M ■ iV)(M 1 / 2 + sdF M ■ N) = M + eN + 0{e 2 ). This implies that 

\\dF M \\ < IIM^H-I^^IIM-^II^ 

Now we apply this bound several times to estimate d p fL and cL/Z . In our situation, 
we have dA = O(l), dB = 0(1). Using M = A and B, we get d(A 1 / 2 ) = 0(\~ 1 / 2 ) 
and d(B 1 / 2 ) = O(l) resp. Using M = B^AB 1 / 2 we get d{B 1 / 2 AB 1 / 2 ) 1 / 2 = 0(A~ 4 / 2 ). 
We now recall that the differential at M of the map M i — > M^ 1 is the linear map 
N i — ► -M^NM- 1 , whose norm is bounded by ||M _1 || 2 . At M = (B 1 / 2 AB 1 / 2 ) 1 / 2 , we 
obtain 

\\d{B 1 ^AB l l 2 )- l l 2 \\ WM^fWdMW =0{\- z/2 ). 

Using M = B 1 / 2 {B 1 / 2 AB 1 / 2 )- 1 / 2 B 1 / 2 we get dL = dM 1 / 2 = 0(A _1 ) and using M = 
^-1/2(^1/2^1/2)1/2^-1/2 we get dL -i = g M i/2 _ o(a~3/4). xhe other estimates are 

straightforward. □ 
Lemma 4.3. The equations of motion in the new coordinates take the form 

x = - ViA(VeI)a; + e 1/2 0(A~ 1/4 <5 + A~ 3 /y ) + 0(e) 
y = V~eHV~eI)y + e^O^H + X'^p 2 ) + 0(e) 
i = 0(V~eS), 

where p = max(||x||, ||y||) is assumed to satisfy p ^ A. The expression for is not useful 
here. 

Proof. The last part of the statement is obvious. We prove the part concerning x, the 
calculations for y are exactly the same. In the original coordinates the vector field ([3]) can 
be written 

s = B(pf){p s - pt(pf)) + 0(\\p s - pl( P f )\\ 2 ) + O(e), 
p s = eA(pf){9 s - et{pf)) + O(e\\0 s - 0*V)f) + 0(eS). 
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As a consequence, we have 



x = L~ 1 B{p s 

+ l- 1 ■ o(||/ - pt || 2 + e ) + e l ' 2 L ■ o(\\e s - e s j 2 + s) 

+ (d p f L~ r ) pf (9 s - 6%) + e-^OVL^V - Pi) 
-L-\d pf et)pf -e-V 2 L(d pf p$)pf. 

We use the estimates of Lemma 14.21 to simplify (recall also that pf = 0(e5)): 

± = L- l B{p s -pl) + e l ' 2 LA(6 s -9t) 
+ (ep + e) + Oie^X-^p 2 + e^X^S) 
+ OiX-hSp 2 ) + 0(A~ 5 / 4 e<5/5) + OiX-hS + A~ 1/4 e 1/2 <5). 



□ 



Lemma 4.4. In i/ie new; coordinate system (x,y,@,I,t), the linearized system is given by 
the matrix 



VeA 0' 

-v/eA 









+ O^X-^r 1 + yfeX'^p + eA- 5 / 4 + ^ 7 ), 



where p = max(||x| 

Proof. Most of the estimates below are based on Lemma 14.21 In the original coordinates, 
the matrix of the linearized system is: 



O(e) 





In our notations we have 



O(e) d 2 Sp sH + O(e) 8 2 pfpS H + O(e) 0" 



O(e) 
0(1) 














O(e) 
0(1) 







+ 0(6e), 





O(e) B- 


hO(e + ^p) 





d 2 pfpS H + O(e) 


0" 




-eA + OieX- 1 /^) 


O(e) 





O(e) 





L = 


O(e) 


0(1) 





0(1) 








































+ 0(Se), 
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In the new coordinates, the matrix is the product 



L 



d(x,y,e,I,t) 
d(0*,p»,9f,pf,t) 



■ L ■ 



d(d s ,p s ,ef, P f,t) 

d(x,y,e,I,t) 



We have 



L/2 


L/2 







0" 









V~ed pfP l + 0(eA- 3 / 4 p) 











7- 1 


































1 



d(e*, P s ,ef,pf,t) 

d(x,y,e,I,t) 
hence 

' d(0 a , P a ,9f,pf,ty 

d(x,y,e,I,t) 

'^BL- l /2 + 0(e\^ 4 ) -^BL- l /2 + 0(eX~ 1 ^) 0(e\~ 3 / A p + e^X' 1 ] 

eAL/2 + OieX-Wp) eAL/2 + OfcA" 1 ^) 

O(v^) O(v^) 







0( 7 - 1 (5e)+ 



e 3/2 0(A -5/4 p + A -l. 






This expression is the result of a tedious, but obvious, computation. Let us just detail the 
computation of the coefficient on the first line, fourth row, which contains an important 
cancellation: 

V~ed% ap sH d pf pt + V~ed 2 pfpS H + 0{eX~^p + e^A" 1 ) 
=V~td pf (d p sH (p*(pf)) +0(eX- 3 / 4 p + eV 2 X- 1 ) = 0(eX^ 3 ^ 4 p + e 3 / 2 A _1 ). 
We now write 



d(x,y,e,I,t) 
[d(es,ps,0f,pf,t) 



and compute that 

"^A + 0(y^A- 3 /V) 
0(eA" 3 / 4 p) 




+0(^A-V 4 7 - 1 ). 



L- 1 







e-V2 L 0(e-V2A-V4) " 

-6-V2L 0(e-V2 A -V4) 

7 

e-V2 

I 



0(^A' 3 /V) 0(eA- 5 / 4 ) 0' 

-^A + 0(^X~ 3 / 4 p) 0(eA~ 5 / 4 ) 

0{y/h/) 0(^7) 
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□ 

In order to prove the existence of a normally hyperbolic invariant strip (for the lifted 
system), we apply Proposition IB.3I to the system in coordinates (x,y,Q,I,t). More pre- 
cisely, with the notations of appendix [Bl we set: 



u = x, s 



y,ci = (@,t),c 2 =1,0 



x n c 



a- - A/2 a+ + A/2 



We fix 7 = yS and a = yj e\/4D, recall that y/eA ^ 2al, by Lemma 14. 1[ We take 
a = Ae~ 1/2 /2, so that 

a_ — A a i + A 



fe V e 

We assume, as in the statement of the Theorem, that < e < eoA 7 / 2 and that ^ 5 < 
y/e X 2 . We can apply Proposition IB. 31 with B u = {u : \\u\\ ^ p} and B s = {s : \\s\\ ^ p} 
provided 

e 1/4 (A- 3 / 4 5 + \- l ' 2 V~e) < P < 2e^ /4 A 5 / 4 . 

It is easy to check under our assumptions on the parameters that such values of p exist. 
These estimates along with Lemma 14.21 imply that 

< 0(4 /4 A) , < 0(4 /4 A 5/4 e 1/2 ). 

Provided that the cylinder C exists, this gives the first set of estimates in Theorem 14.11 
Let us check the isolating block condition. By Lemma 14.31 we have 

x ■ x > 2a\\x\\ 2 - \\x\\ Oie^X-^H + e^A^V + e) 

if x 6 B u ,y € B s . If in addition ||x|| = p, then 



A- 3 / 4 ^ e J /4 | 



a-3/v 2 < 2 £ y 4 i 



1/4, 



hence 



x ■ x ^ 2a||x|| 2 — ||x|| 2 ej!/ 4 0(\/eA) ^ a|M| 2 

provided eo is small enough. Similarly, y ■ y — a||y|| 2 on B u x dB s provided eo is small 
enough. Concerning the linearized system, we have 

L uu = yfeA + 0( v ^5A- 1 /4 7 -i + yfe\-^p + e A~ 5 / 4 + yfa) 

= yfeA + 0(el /A V7\) ^ al, 

L ss = -yfeA + 0(eo /4 \/eA) < -al 

on B u x B s x O r . These inequalities holds when eo is small enough because y/eA 2al 
and \/eA ^ 0(a). Finally, still with the notations of Proposition IB. 31 we take 



m 



0{y/lSX 



-V4 7 -l 



+ yfe\- z ' A p + eA" 5 / 4 + yfe 1 + \fe8/a) 



1/4, 
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If eo is small enough, we have 4m < a hence 

K < 2m/a < O^ 4 ) < 2" 1 / 2 , 

and Proposition IB.3I applies. The invariant strip obtained from the proof of Proposition 
IB. 31 does not depend on the choice of p. It contains all the full orbits contained in 



{x : |N| < 4 /4 A" 5 / 4 } x {y : ||y|| < e^A" 5 / 4 } 



x K x 



a_ - A/2 a + + A/2 



hence all the full orbits contained in V, as defined in the statement of Theorem 14.11 The 
possibility of taking p = e 1 ^ 4 (A _3 / 4 5 + \~ l l 2 ^fe) now implies that the cylinder is actually 
contained in the domain where 

IMIJvK^a^ + a-vv?). 

Moreover, with this choice of p and using that K = 0(m/\/eA), we can obtain an improved 
estimate of the Lipschitz constant K: 

K = 0{V~5\- 3 ^ + P A- 5 / 4 + V~^- 7/4 ) 

= 0(Vd\- 3 /* + e- 1/A 5\~ 2 + e" 1/4 v^A" 7 / 4 + ^A~ 7 / 4 ) 
= 0(^A- 3 / 4 + V6X- 1 + e- 1/4 ViA- 7 / 4 ) 
= O(V^A- 1 + e - 1/4 ViA- 7 / 4 ). 

Observe finally that, since the system is l/7-periodic in © and 1-periodic in t, so is the 
invariant strip that we obtain, as follows from Proposition IB.2I We have obtained the 
existence of a C 1 map 



w c = (w c u ,w c s ) : (9,J,i)GRx 



a_ - A/2 a + + A/2 



t n-l x R n-1 



which is 2-ftT-Lipschitz, l/7-periodic in and 1-periodic in t, and the graph of which is 
weakly invariant. 

Our last task is to return to the original coordinates by setting 

e s (9f,pf,t) = 9t(pf) + h(pf) ■ « + 0(7f / ,e- 1/ V,t) 

All the estimates stated in Theorem 14.11 follow directly from these expressions, and from 
the fact that ||d?i; c || ^ 2K. This concludes the proof of Theorem 14. 11 □ 
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5 Localization and Mather's projected graph theorem 

We now study the system in normal form N t = Hq + eZ + eR from the point of view of 
Mather theory. We study the normal form system N e = Hq + cZ + ei? on the neighborhood 
of the set {p = p*(p^),p^ G [a_,a+] C [a m j„, a max ]}. We assume that Z satisfies the 
generic conditions [G0]-[G2] and that ||-R||c 2 ^ ^- Recall that there exists a partition of 
[a m i n ,a max ] = Uj=i[°ji a j+i]> such that forp-f G [aj — A, a J+ i + A] the function Z(9 s ,p s ,pf) 
as a nondegenerate local maximum at 0|. It is clear that we can restrict this partition to 

[a_,a+]. We abuse notation and still write [o_,a + ] = IJj=i[ a i' a j+i\- 

We first point out the following consequences of the genericity conditions [G0]-[G2]: 
there exists < b < A/4 depending on Hi such that 

[Gil 

Z(9°(pf),p*(pf)) - Z{6 s } p*(pf)) > b\\9 s - (9?(p/)||, 
for each p$ G [oy + 6, Oj+i — 6]. 
[G2'] For pf G [oj+i — b, aj+i + 6], j = 0, • • • , s — 2, we have 

max{Z(^,p*(/)),Z(^ +1 ,p*(/))} - Z{9 s Mpf)) 

>bmin{\\e s -e%\\e s -e s j+ i\\y. 

In the first case, the function Z has a single non-degenerate maximum, which we will 
call the "single peak" case, while the second case will be called the "double peak" case. 
The shape of the function Z allows us to localize the Aubry set and Mane set of the 
Hamiltonian N e . 

According to Theorem 14.11 for each [aj — A/2, a J+ i + A/2] there exists 
Xj = {(B s ,p s ) = (e*,P;)(^,/,i)); pf G [a, - ^,a i+1 + ±], (0 f ,t) G T x T}, 
which are maximally invariant set on Nj := {(9,p,t);pf G [aj — ^,a J+ i + £], \\(9 s ,p s ) — 

These information allows us to study the Mather set, Aubry set and Mane set of the 
Hamiltonian N e . 

Theorem 5.1 (Localization). For N e = Hq + eZ + eR such that Z satisfies [G0]-[G2], 
then there exists e$, So and < pi < pi such that for < e < eo and < 5 < 5o the 
following hold. 

1. For any c = (p*(c-' ), & ) such that cf G [aj + b, Oj+\ — b], Af(c) is contained in 
{(0,p,t),\\ P -c\\^6AV^,\\9 s -e s j (pf)\\^p 2 }. 
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2. For c = (pt{cf),cf) such that c$ G [a-j+i — b,aj + ± + b], we have that A~nX c ) ^ s 
contained in 

{(d,p,t), \\p - c\\ < 6A^min{\\9 s - 9^)1 \\9 S - 9 s j+1 (pf)\\} ^ p 2 }. 

Apply the statements of the previous theorem with Theorem 14, 1\ we may further 
localize these sets on the normally hyperbolic cylinders. Moreover, locally these sets are 
graphs over the 6$ component, which is a version of Mather's projected graph theorem. 

Theorem 5.2 (Mather's projected graph theorem). For any N e such that Z satisfies 
[G0]-[G2j, there exists 5o and eo depending on b, n and r such that for 5 ^ 5q and e < eq 
we have: 

1. There exists < p 2 < pi such that for for c = (p^ ), cf ) with c* G (a-,- + b, aj + i — b) 
the Marie set N c is contained in the normally hyperbolic cylinder Xj . 

Moreover, let ir 9 f be the projection to the 9? component, we have that -Kqj\A c is 
one-to-one and the inverse is Lipshitz. 

2. For a G [aj+i — b, aj+i + b], we have that A c C Xj U Xj + \. 

TTgf\A c r] Xj and ir 9 f\A c r) Xj + i are both one-to-one and have Lipshitz inverses. 

For the rest of this section, we will derive various estimates of quantities and set arising 
from Mather theory. We deduce Theorem 15.11 and Theorem 15.21 from these estimates. 

5.1 Vertical estimates 

We derive estimates on the Mather sets of a general Hamiltonian H(t,9,p), depending on 
e, under the assumptions that 

I/A ^ d P pH ^ AI 
in the sense of quadratic forms, and 

\\d e H\\ c i s£ 2e. 

Note that both Hamiltonians H t and N e satisfy these assumptions. The main result in 
this section is: 

Proposition 5.3. We assume that e ^ 1. For each cohomology c G W 1 and each Weak 
KAM solution u of H t at cohomology c, the set Z(u, c) is contained in a 36 A^fi- Lipshitz 
graph, and in the domain \\p — c\\ ^ QAyJne. 

It is useful to use the Lagrangian L{t,9,v) associated to H. Recalling the expressions 
8 vv L(t,9,v) = (d pp H(t,9,d v L(t,9,v)y\ 
d Bv L(t,9,v) = -d dp H(t,9,d v L(t,9,v))d m L(t,9,v) 
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and 

d e0 L(t, 0, v) = -dg 9 H(t, 9, 8 v L{t, 0, «)) - d 9p H(t, 0, d v L{t, 0, v))d dv L(t, 6, v), 
we obtain the estimates 

ll^wLHcro < A \\devL\\c ^ 2Ae , \\deeL\\c° ^ 3e 

when e < eo(A). 

We recall the concept of semi-concave function on T". A function u : T n — > R is 
called K-semi-concave if the function 

x i — ► u(x) - K\\xf/2 

is concave on R n , where u is seen as a periodic function on R n . It is equivalent to require 
that, for each £T n , there exists a linear form I on R n such that the inequality 

u(0 + y)^u(0) + l-y + K\\y\\ 2 /2 

holds for each y G R n . The following Lemma is a simple case of Lemma A. 10 in [Belj : 
Lemma 5.1. If u : T n — > R is K -semi-concave, then it is (K \/n)-Lipshitz. 
Proof. For each i£T", there exists l x E R n such that 

u(x + y) ^ u{x) + l x -y + K\\y\\ 2 

for all y G R n . By applying this inequality with y = (±1,0,0, ••• ,0), we conclude that 
the first component (^3;)i of l x satisfies |(^)i| ^ K. Similar estimates hold for the other 
components of l x , and we conclude that that \\l x \\ ^ Ky/n for each x, and thus that u is 
-KT-^/n-Lipshitz. □ 
We will need the following regularity result of Fathi: 

Lemma 5.2. Let u and v be K-semiconcave functions, and let I C T n be the set of points 
where the sum u + v is minimal. Then the functions u and v are differentiable at each 
point of I, and the differential x i — > du(x) is 6K -Lipshitz on I. 

Let us recall that the Weak KAM solutions of cohomology c are defined as fixed points 
of the operator % : C(T n ) — > C(T n ) defined by 

%(u)(0):=minu(7(O))+ / L(t,j(t),j(t)) + c ■ j(t)dt, 

7 JO 

where the minimum is taken on the set of C 1 curves 7 : [0, 1] — > T n satisfying the final 
condition "f(T) = 0. 

Proposition 5.4. For each c G R n , each Weak KAM solution u of L + c-v is 6A^/e-semi- 
concave and 6A^/ne-Lipshitz. 
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Proof. For each T G N and 9 €T n , we have 

u(fl) = minu( 7 (0)) + / L(t,7(t),7(t))+c-7(t)dt, 
t Jo 

where the minimum is taken on the set of C 1 curves 7 : [0, T] — > T n satisfying the 
final condition j(T) = 9. Let @(t) be an optimal curve in that expression, meaning that 
G(t) = 9 and that 

u(9) = u(9(0)) + [ L(t, 9(t), @(t))dt. 
Jo 

We lift (and the point 9 = O(T)) to a curve in M n without changing its name, and 
consider, for each x E R n , the curve 

O x (t) ■= @(t) + tx/T, 

so that Q X (T) = 9 + x. We have the inequality 

u(9 + x)-u(9)^ [ L(t,@ x (t),@ x (t))-L(t,9(t),0(t)) + c-x/Tdt. 
Jo 

The integrand can be estimated as follows: 

L(t,e x (t),e x (t)) ^ L(t,e(t),e(t)) 

+ dgL(t, @(t),@{t)) ■ tx/T + 8 v L(t, @(t),@(t)) ■ x/T (5) 
+ 3e\tx/T\ 2 + 2Aet\x/T\ 2 + A\x/T\ 2 /2 

Integrating, and using the Euler-Lagrange equation, we conclude that 

u{9 + x)- u(9) ^ (c + d v L(T, G(T), 9(T)) • x + (eT + e + l/2T)A\x\ 2 

for each T G N. Taking T £ [1/2-^/e, l/\/e] (this interval contains an integer since e ^ 1), 
we obtain 

u(9 + x)- u{9) ^ (c + 9„L(T, 9(T), 9(T)) • x + 3A^\x\ 2 . 

This ends the proof of the semi-concavity. The Lipshitz constant can then be obtained 
from Lemma 15. II □ 
Let u be a weak KAM solution, and let u be the conjugated dual weak KAM solution. 
Then the set Z(k, c) can be characterized as follows: Its projection X(n, c) on T n is the set 
where u = u, and 

I(u, c) = {(x, c + du(x)),x G X(u, c)}. 

Since — u is semi-concave, it is a consequence of Lemma 15.21 that the differential du{x) 
exists for x G X(n, c). Moreover, we can prove exactly as in Proposition 15.41 that —u 
is 6^4-^/e-semi-concave. Lemma 15.21 then implies that the map x 1 — > du(x) is 36A^/e- 
Lipschitz on X(n, c). Moreover, du(x) is bounded by the Lipschitz constant &Ay/ne of u. 
This ends the proof of Proposition 15.31 □ 
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5.2 Horizontal localization 

For 6 s <E T n_1 and r > 0, let D(9,r) denote the closed Euclidean ball centered at 9 s with 
radius r. 

Proposition 5.5. Let N e = Hq + eZ + eR with Z satisfying [GO] — [G2], then there exists 
k > depending only on A, b and n, eo > depending only on A and 5 such that for each 
e e]0, eo[ and c G T(e) we have the following results on the projected Mane set and the 
Aubry set. 

1. If c = (p%{cf),cf) with c* G [aj +b,aj + i — b], then for each weak KAM solution u of 
N e at cohomology c we have that 

l(u,c) C D (0t(cf),K6ty x T, 

as a consequence, we have 

M Ne (c) C D (0 9 j{c*),K5*} x T. 

2. If c = ip%(c*),c*) with c* G [oj+i — Oj+i + b], then 

AnM C (l> (^(</),/eJ2) xl)u(i) (^+i(c / ),k<5s) x T) . 

The Lagrangian N*(t,9,v) associated to N e will play a central role in the proof. We 
write it 

N*{t, 6, v) = L {v) + eZ{9 s , d v L {v)) + eL 2 (t, 9, v, e), 
where Lq is the Legendre dual of Hq. We have 

21 1 A ^ d vv L ^ AI/2 

in the sense of quadratic forms. 
Lemma 5.3. We have the estimate 

miH 2 s$ L 2 (t,9,v) < Ae + supH 2 . 
Proof. Let us first consider the truncated Hamiltonian 

H(t,9,p) = H (p) + eZ(9 s ,p) 
and the associated Lagrangian L(9 s ,v). We claim that 

L {v) - eZ(9 s , dL (v)) < L(9 S , v) < L (v) - eZ(9 s , dL (v)) + (6) 
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In order to prove the left inequality, we write 

L(9 s ,v)= S np[p-v-H (p)-eZ(9 s ,p)} 
v 

> dL (v) • v - H (dL (v)) - eZ(9 s , dL (v)) 
>L (v)-eZ(6 s ,dL (v)). 

The right inequality follows from the following computation: 

L(9 s ,v)=snp[p-v-H (p)-eZ(9 s ,p)} 

p 

p-v- H (dL (v)) -v(p- 8L (v)) - \\p - dL (v)\\ 2 /A 

eZ(e s ,dL (v)) + e\\p-dL (v) 
v ■ dL (v) - H (dL (v)) - eZ(9 s , dL (v)) 

+ e\\p-dLo(v)\\-\\p-dL (v)\\ 2 /A 



^ sup 

p 



^ sup 

p 



< L (v) - eZ(0 s ,dL o (v)) + sup[ey - y 2 /A)} 

«S L (v)-eZ(9 s ,dL (v)) + Ae 2 . 
Now we have estimated L, we observe that 

H(t,9,p)-esupH 2 «S N e (t,9,p) ^ H(t,9,p) - einf H 2 

from which follows that 

L(t, 9,v) + e inf H 2 ^ L(t, 9, v) ^ L(t, 9,v) + e sup H 2 , 

which implies the desired estimates in view of ([6]). 
Let us now estimate the a function of N e : 



□ 



Proposition 5.6. The a function of Mather is estimated at the points c S T in the 
following way: 

H (c) + eZ(e s (c), c) - e max L 2 (t, 9, dH (c)) < a(c) < 

^H (c) + eZ(9 s (c),c)-e min H 2 (t,0,c) 

(t,6»)eT n + 1 

thus 



H Q (c) + eZ(9 s (c),c) - e\\H 2 \\ c o - Ae 2 < a(c) < H (c) + eZ(9 s (c),c) + e\\H 2 \ 



c» 
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Proof. We have 

a(c) < max Hit, 9, c) ^ Hn(c) + e max Z(6 S , c) — e min H 2 it,0,c) 

(t,d) (t,e)eJ"+ 1 

which is the desired right hand side. On the other hand, let us set oj = dH$(c) € W l 
and observe that c = 8Lq{u). We can consider the Haar measure fx of the torus T x T x 
{©•^(c)} x {oj}, this measure is not necessarily invariant but it is closed. We thus have 

a(c) ^ c • oj — J N*dfi = c ■ oj — Lq{oj) + eZ(Q^(c), c) — e J L 2 d[i 

^ H (c) + eZ(O f (c),c) - e max L 2 (t,6,u) 

(t,e)eT«+ 1 

□ 

Lemma 5.4. For eac/i c £ T, have the estimates 

N*{t,9,v)-c-v + a{c) ^ \\v - dH (c)\\ 2 /(2A) - eZ c (6 s ) - ct] (7) 
iV*(t, 0, v) - c ■ v + a{c) ^ A\\v - dH (c)\\ 2 /2 - eZ c (9 s ) + erj (8) 

where Z C (9 S ) = Z{0 S , c) - max es Z(9 S , c) and 

r] = 2\\H 2 \\ c o + ( 2A + A 3 )e. 

Proof. It is a direct computation: 

N*{t,9,v)-c-v + a{c) ^ 

s: L (v) -c-v + H (c) - eZ(6 s , d v L (v)) 

+ emaxZf^c) + Ae 2 + esupH 2 - emmH 2 (t,9,c) 

e s 

^ A\\v - dH (c)\\ 2 /A - eZ c {9 s ) + e\\dL (v) - c\\ + Ae 2 + 2e|| J ff 2 ||c° 
^ A\\v - dH (c)\\ 2 /A - eZ c {9 s ) + Ae\\v - dH (c)\\ + Ae 2 + 2e||F 2 || c o 
^ A\\v - dH (c)\\ 2 /2 - eZ c {9 s ) + 2Ae 2 + 2e\\H 2 \\ c » 
N*(t,9,v)-c-v + a(c) ^ L (v) - c ■ v + H (c) - eZ(8 s ,d v L (v)) 
+ maxZ(9 s ,c) -2e\\H 2 \\ c o - Ae 2 

6 s 

> \\v - dH (c)\\ 2 /A - eZ c (9 s ) - e\\dL (v) - c\\ - 2e\\H 2 \\ c o - Ae 2 
^ \\v - dH (c)\\ 2 /A - eZ c {9 s ) - Ae\\v - dH (c)\\ - 2e\\H 2 \\ c o - Ae 2 

> \\v - dH (c)\\ 2 /(2A) - eZ c (9 s ) - Ae 2 - A\ 2 - 2e\\H 2 \\ c0 



□ 



We can now estimate the oscillation of a weak KAM solution near 9? and 0j +1 . 
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Lemma 5.5. Let u(t, 9) be a weak KAM solution at cohomology c. 

1. If c = (pl(cf),cf) with G [cij + b,aj + i — b], then for any fa, 9i), fa, 9 2 ) € T x 
D(e s j(cf)) x T 

u(t 2 , 6» 2 ) - ufa, 6»i) < 4n VnAe 

where r\ = ^Arj/b. 

2. If c = {vl{c f ),cf) with of G [aj+i — 6, cij+i + 6], i/jen /or either (ti,9i), (t 2 ,9 2 ) G 
T x L»(^(cO) x T or fa, fa, #2) G T x D{9 s j+1 { C f)) xT, 



u(t 2 ,9 2 ) - u(h,9i) < AnVnAe. 
Proof. Using ||^||c2 ^ 1, we have that 

4(fl*)>-^r-^)ii 2 . 

We take two points fa, i = 1 or 2 in this domain, and consider the curve 

9( t) =9 1+{ t- ^hzh+m^*^^ 

t 2 -h+T 

where T G N is a parameter to be fixed later, and where U G [0, 1[ and 9{ G [0,1 [ n are 
representatives of the angular variables U,0i, and [u>] G IT 1 is the component-wise integral 
part of u. Note that 9(h) = 9\ and 9(t 2 +T) = 9 2 , hence 



ft-2+T 

u(t 2 ,9 2 ) -ufa,0i) ^ / L{t,9(t),9(t))-c-9(t) + a(c)dt 
Jti 

rh+T 

/ A\\9-dH {c)\\ 2 /2 + eZ c {9 s {t)) + e V dt 
f i2+T 2An 2/ 

r*2+T 2An 2 , 



This inequality holds for all T G N, in particular, we can choose T G N so that 

IfiA m ~ 

< T + t 2 - h < 2 





and obtain 

u(t 2 , 9 2 ) - ufa, 61) < 4riVnAe. 

□ 
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Up to now, we used that ||Z||c2 ^ 1, but we used no information on the shape of Z. 
Now we use properties [Gl'] and [G2'] to prove Proposition 15.51 



5.2.1 The single peak case 

This concerns the first statement of Proposition 15.51 where 

c = p*(t/), G G [dj + b, a,j+i — b]. 
By [Gl'] the function Z(6 S , c) as a single peak at 9?, as a consequence 

Let 9{t) : M. — > T m be a curve calibrated by u. Then the function 

1 i— ► L(i, 9{t),0(t)) - c ■ 9(t) + a(c) 

is integrable. Since 

L(i, 6, v) > -eZ c (# s ) - e?7 ^ e&r 2 - er? > e&r 2 /4 

if # s does not belong to D(0j,n), we conclude that the set of times t for which 9 s (t) does 
not belong to D c {r{) has finite measure, and is an open set. Let ]ii,i2[ be a connected 
component of this open set of times. Then 9 s (ti) and 9 s {t2) belong to D(9j,ri) hence 

ti 

L(t, 6(t),d(t)) - c • 9{t) + a{c)dt = u(t 2 , 9(t 2 )) - u(h, 9(t x )) < AnVnAe. 

h 

Now let ro be the maximum of the distance ||# s (i) — G s (c)||, assume that ro ^ 2r\. Let t$ 
be the smallest solution of the equation — © s (c) || = ro in ]*i , *2 [, and let < t\ be 

the greatest solution of the equation ||# s (i) — O s (c)|| = tq/2 in [i i , ^4] . Note that 

I * L{t, 9(t), 9(t)) - c ■ 9(t) + a(c)dt ^ AnVnAe 
Jt 3 

because the integrand is positive on ]ii , ^2 [■ We conclude that 
Jt 3 

and, using the Cauchy- Schwartz inequality, that 

fti 1 / /"*4 \ 2 chr 2 



\9 s (t)\\ 2 /(2A) + ebrl - erjdt < AnVnAl, 



Jf r(t)|| 2 /(2^) + ebrl - e V dt > ^ _ ^ (jf 11^(011*] + - *a) 



r 2 , e&r 2 , , , rl^/be 

" 2A(t 4 - is) 2 1 3; " 2^/Z 



(9) 
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Finally, we obtain 

r^y/be 
2\/A 

or equivalently 



^ 4r±V nAe. 



2 , i— 8AV2n _ 

?~o ^ T\aA\/ n/b = — - — y/r). 

5.2.2 Double peak case 

We now turn to the case of 

c = p*(c / ), d G [cij+i - b,ctj + i + b], 
where the function Z(8 s ,c) has two potential maxima. It follows from [G2'] that 

z c {e s ) < -b (mm{\\e s - \\e° - e° +1 \\}) 2 . 

Let 8 = (0gX) G A Ne (c) be where the function (of 8) min{||0 s -0|||, ||6> s -(9j +1 ||} achieves 
its maximum. This is possible since AnX c ) ls a compact set. Since /i(#o,$o) = 0, there 
exists an increasing sequence of integers n k and absolutely continuous curves 8 k ■ M — > T n 
satisfying 8 k (0) = 9$ and 8 k (t + n k ) = j k (t), and 

lim / L(t,8 k (t),8 k (t))-c-8 k (t) + a(c)dt = 0. 

k — >oo J Q 

Similar to the first case, L(t,9,v) ^ ebrf/4 for 9 s D(8 S j,r\) U D(8j +1 ,ri), we conclude 
that for sufficiently large k, 8 k (M) must intersect D(8j,ri) U D(0| +1 ,ri). Without loss of 
generality, we may assume that it intersect D(9j,rx). 

Let h = min{ti < 0,6 k (ti) G D{9 s j ,r 1 )} and t 2 = min{0 < t 2 ,9 k {t 2 ) G D(0j,ri)}. We 
first study the action of 9 k on the interval [0, t 2 ]. If 9 k ([0, t 2 ]) does not intersect D(8j +1 ,ri), 
write *3 = £4 = t 2 , otherwise, Write £3 = min{0 < t 3 ^ £2j7(*a) £ D(9j +1 , n)} and 
t 4 = max{t 3 < t 4 < *2,7(t 4 ) G L>(0j +lJ ri)}. 

We still use ro to denote the maximal distance min{||$g — ll^o — ^j+ill}> assume 
that ro ^ 26i. Let £5 G [0,i3] be the smallest solution in [0, £3] such that d(8 s {t§)) = ro/2, 
then by the same calculation as in ([9|), 



t5 L{t, 8 k (t),8 k (t)) - c -8 k (t) + a(c)dt > 
2\IA 

Furthermore for any weak KAM solution u(8,t) 

rt4 



L-c- 8 k {t) + a{c)dt > u(U, 8 k (U)) - u(t 3 , 8 k (t 3 )) ^ -AnVnAe, 



*3 
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while the integrand is nonnegative on both [^,£3] and [£4, £2]- We conclude that 

rte r 2 rz 

/ L - c ■ 9 k (t) + a(c)dt ^ -AriVnAe. 
Jo 2VA 

The same estimate can be made for the action on the interval [ti,0]. In addition 

rn k +t! _ 

/ L- c - 9 k (t) + a(c)dt ^ u(n k + h,0 k (rik + h)) - u{t2,9 k (t 2 )) ^ -ArWnAe, 
Jt 2 

note that t 2 < n k + 1\ and that k is n k periodic. 
Finally we conclude that 

rrik+ti r 2 \fbe 

/ L — c - 9 k + a(c)dt ^ -^-7= — \2bivnAe. 
J tl 2VA 

Let k — > 00, the integral on the left hand side approaches 0. We obtain 

r^Vbe 



2y/A 

and 



^ 126i VnAe 



2 2AA^fn 
r ^ - h ^7. 

We choose eo sufficiently small such that rj ^ 2[|Ho||c , choose k = \J ' ^AV2n an( ^ ver jfy 
that we have proved the statements of Proposition 15,51 in both cases. □ 

Before moving on we point out that the estimates in the double peak case indeed implies 
that the curves 6 k , which are not calibrated, can be localized in the limit of k — > 00. We 
state it in the following lemma for future use. 

Lemma 5.6. For the double peak case, i.e. c = p*(cf) with c* € [a^+i — b,aj + \ + b], 
consider any 9q G An c (c), ^ n k be an increasing sequence of integers, 6 k = (0 k ,Q k ) '■ 
M — > M be a sequence of n k — periodic absolute continuous curves such that 7(0) = 9q and 



f n k 

lim / L(t,9 k ,9 k ) - c- 8 k + a(c)dt = 0, 



Am I 

tc- 
theri there exists K S N such that for all k > K 

max min{||0|(i) - 9% \\e s k {t) - < 2k&\. 

Proof. Fix a curve fc , write d{6 k (t)) = min{ ||6>|(t)— 6>| || , \\9 s k (t)-9 s j+l \\}. Let r be where the 
maximum of d{6 k {t)) is reached. Consider the shifted curve 9' k (t) = 9 k (t—r), the arguments 
in section 15.2.21 go through with 6 k replaced with 9' k and ro replaced by maxd(9 k (t)). We 
have that 

lim max d(9 k (t)) ^ k6* 

k — >oo tgR 

and the lemma follows. □ 
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5.3 Proof of Theorem I5TT1 



Proposition l5.3l provides the vertical part of the localization, while Proposition 15 , 51 provides 
the horizontal localization, with p 2 = k<5* . Clearly we can choose <5o small enough such 
that p%< p\. □ 

5.4 Proof of Theorem 15.21 

For the first case, where c = p*(cf) with & G [aj + b,a,j+i — b]. For a sufficiently small 
choice of eo, Theorem 15.11 implies that A/j\r e ( c ) C Vj, where Vj = {{9,p,t);p* G [a,j — 
b,aj + i + b], \\(9 s ,p s ) — (9j,pt)\\ ^ Pi} was defined in Theorem 14.11 Since Vj is maximally 
invariant and that Mn € (c) is an invariant set, we conclude that A/jv £ (c) C Xj. 

For the second case, where c = p*(cf) with c-* G [oj+i — fr, Oj+i + 6], we can similarly 
claim that AnX c ) C V,- U Vj+i, moreover AnX c ) ^ Vj an d Av e (c) H Vj+i must both be 
invariant, and hence ^4at £ (c) n Vj C Xj and An c (c) n Vj+i C Xj + \. 

In order to prove the projection part of Theorem 15.21 let us consider a Weak KAM 
solution u of N e at cohomology c. Let (U,9i,pi),i = 1,2 be two points in X{u, c). We shall 
denote by the same symbol k various different constants which are independent of 5 and 
e. By Proposition 15.31 we have 

\\P2-Pi\\ ^ 36AV~e\\9 2 - e x \\ < 36^(||^ - o{\\ + \\9 S 2 - 9{\\). 
Assume that these two points belong to one of the NHICs Xj, we also have 
P s 2 ~ 0f II < (K$/V~e)(\\9 f 2 - 9{\\ + \\ P2 - Pl ||). 

We get 

\\p 2 - Pl \\^ K V-e\\9 f 2 -9{\\ +K 5\\p 2 -p 1 \\ 
thus, if 5 is small enough, 

\\P2-PiH">ft\\4-rt\\ 

and then 

m-ei\\<w,/e)\\0i-e{\\. 

We have proved that the restriction to T(u, c) of the coordinate map 9? has a Lipschitz 
inverse. 

Note that the Mane set A/jv e (c), as well as the components of the Aubry set ,A/v £ (c)n Vj 
and -4at e (c) n Vj+i are both contained in some I(u, c), since we have just proved that they 
belong to NHIC, their projection to the Q* component has a Lipshitz inverse. □ 

6 Variational Construction 

More detailed information on these sets can be obtained, if we are allowed to make an 
additional perturbation to avoid degenerate situations. 
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Theorem 6.1. Let N e = H$ + eZ + eR be such that Z satisfy the genericity conditions 
[G0]-[G2] and that the parameters e and 5 is such that Theorem \4-l\ applies. Then there 
exists arbitrarily small C r perturbation eR' of eR, such that the following hold for N' e = 
H + eZ + eR': 

1. There exists a partition of [a_,a + ] into Ui=o[%>^j+lL which is a refinement of the 
partition {[aj, aj+i]}. Each [aj, still corresponds to an invariant cylinder Xj. 
We have that for a G (a,j,a~j+i), the Aubry set An^(c) is contained in Xj; for 
& = cij+i, An^(c) has nonempty component in both Xj and Xj + \, if Xj ^ Xj + \. 

2. The sets An^(c) n Xj, when nonempty, contains a unique minimal invariant proba- 
bility measure. In particular, this implies that An{, (c) = A/jv' ( c ) f or c ^ 7^ ®j f or an V 
3 ■ 

3. An immediate consequence of part (2) is the following dichotomy, for cf ^ ~dj, j = 
1, • • • ,1, one of the two holds. 

(a) A c = N c and 7TgfA c = T. In this case, A c is an invariant circle. 

(b) TT fJ\f c C T. 

Using the information obtained from the normal form system A^^, we now return to 
the original Hamiltonian H t . Using the symplectic invariance of the Mather, Aubry and 
Mane set developed in [Be2j , we have that the same conclusion as in Theorem 15.21 and 
Theorem 16,11 can be drawn about H e . 

Theorem 6.2. Let H e = Hq + eH\ such that the resonant component of Hi satisfy con- 
ditions [G0]-[G2]. There exists eo > and an interval [a_,a+] C [a m j n , a max ] depending 
only on Hq and H\, and for each < e < eo there exists arbitrarily small C r perturba- 
tion H' e of H e , such that the conclusions of Theorem \5.2\ and Theorem I6.il holds for the 
Hamiltonian H' e at cohomologies c = p*(c^), where € [a_,a+]. 

These information on the Mane set allow to use the variational mechanism of [BelJ. 
Let us denote by T(e) the set of cohomology classes c G T such that G [a_,a+]. We 
would like to prove that each cohomology c G T(e) is in the interior of its forcing class 
in the terminology of [Bel] , which implies that all the cohomology classes in T(e) are 
contained in a single forcing class. By proposition 5.3 in [Bel] , we could conclude the 
existence of an orbit (9(t),p(t)) of H € such that p(0) = c and p(T) = d for some T G N. 
Note that this implies the existence of various more complicated orbits, see [BelJ. 

In order to carry out this program, we denote by To(e) the set of cohomology classes 
c G r(e) such that the set 0f(I(c, u)) is properly contained in T for each weak KAM 
solution u at level c. By Theorem 0.11 in [Bel] , each cohomology c G ro(e) is in the 
interior of its forcing class. 

Let T2(e) denote that set of c G T(e) such that the Aubry set A{c) has exactly two 
static classes. In this case the Mane set M(c) D A(c). To ensure that c G ^(e) is in the 
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interior of its forcing class, some further degeneracy conditions are needed. To be more 
specific, let T^c) denote the set of c G ^(e) such that the set 



M{c) - A{c) 

is totally disconnected. This can also be stated in terms of barrier function. Let 9q and 
9\ be contained in each of the two static classes of A(c), we define 

bt(9) = h c (0 ,6) + h c {0,0 1 ) 

and 

b-{0) = h c {9 1 ,0) + h c (0,9 ), 

where h c is the Peierls barrier for cohomology class c. Then T^e) is the set of c G ^(e) 
such that the minima of each b + and b~ outside of A(c) are totally isolated. 

We call Ti (e) the set of cohomology classes c G T such that there exists only one weak 
KAM solution u at level c, and 6? (l(c,u)) = T. Note then that 

M{c) = A(c) = I(c, u) 

is an invariant circle. We have To(e) n ri(e) = for each e g]0, eo[, by definition. We first 
consider the covering 



e = {ef, eie s 2 ,---, e s n _ x ) —>• m = (o f , 20?, 2 v ■ ■ , 9 s n _ ± ). 

This covering lifts to a a symplectic covering 

(9,p) = (9,pf,plp S 2 , . . . ,K-l) — > £(M = (m,P f ,P S l/2,P S 2, ■ ■ ■ ,Pn-l), 

and we define the Lifted Hamiltonian H = H o S. It is known that 

A k {c) = ~- 1 (A k {c)) 

where c = £*c = (c/, cf /2, c|, . . . , c*_ 1 ). On the other hand, the inclusion 

N fl (c)DZ- 1 (N fl (c))=Z~ 1 (A il (c)) 

is not an equality for c G Ti(e). More precisely, for c G Ti(e), the set ^t^-(c) is the union of 
two circles, while M k {c) contains heteroclinic connections between these circles. Similarly 
to the case of ^(e), we call r*(e) the set of cohomologies c G Ti(e) such that the set 
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is totally disconnected. Alternatively, we can chose a point 9q in the projected Aubry set 
A(u, c) of H, and consider its two preimages 9q and 9\ under ^. We define 

b+(e) = h(9 ,6) + h(e,e 1 ) 

and 

b-(p) = h0x,o) + h(e,o o ) 

where h is the Peierl's barrier associated to H. r*(e) is then the set of cohomologies 
c G Ti(e) such that the minima of each of the functions b^r located outside of the Aubry 
set Afj(c) are isolated. 

The following theorem is proved in [BclJ. 

Theorem 6.3. If c and d belong to the same connected component o/To(e)ur|(e)ur2(e), 
then there exists an orbit (9(t),p(t)) and of H e a time T G N such that p(0) = c and 
p(T) = d. 

We have proved the main result provided T(e) = To(e) U r|(e) U T*(e). 

Theorem 6.4. Let H e be a Hamiltonian such that Theorem \6.S\ holds, then there exists 
an arbitrarily small C r perturbation H" to H' e , such that for the Hamiltonian H' e ' we have 
that T(e) = r (e) U r|(e) U I^(e). 

We note that the conclusions of Theorem 16.11 already implies that {c G T(e),c* ^ 
aj,j = 2,---l- 1} C r (e) UTi(e), while {c G T(e),tJ = a,,j = 2, • • • I - 1} C r 2 (e). 
In other words, T(e) = To(e) U r x (e) U T2(e). It suffices to prove that Ti(e) = T\(e) and 

r 2 (e) = rs(e). 

For the rest of this section, we prove Theorem 16.11 and Theorem 16.41 
6.1 Local extension of a(c) and AnX c ) 

Consider the normal form system N e and pick c = p*(cf) with G [a_,a+]. For such a 
c the function Z(9 s ,c) has a single peak. It follows from Theorem 15.21 that the Aubry set 
Ai\r e (c) (which is a subset of J\Fn c (c)) is contained in a single NHIC Xj and the projected 
graph theorem holds. For the rest of the cohomology classes, the double peak case, the 
picture is less clear as An c (c) are contained in the union of two NHICs. To get a more 
precise picture, we will locally extend the set function (of of) 

•An c ((pt(c f ),C f ))\ [a . +b ,a J+1 -b] 

from [a,j + 6,07+1 — 6] to [aj — 2b, Oj+i + 2b]. The extended local Aubry set will still 
be contained in Xj. These definitions are inspired by Mather's definitions of relative 
a— function and Aubry set. 
Let 

P0 = min ||0j(/)-0J +1 (p/)||/3. 

pJ e[aj + i — 2o,aj + i+26] 
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It follows from properties [Gl'] and [G2'] for pf G [aj — 26, aj+i + 26], 

z(ei, P i,pf) - z(e s ,pt, P f) > b\\e s - 9 S \\ 2 

for \\9 S — 6j\\ ^ po- By choosing a smaller 5 if necessary, we may make sure po > p\, 
where pi was defined in Theorem 14.11 We write Uj(p) = {\\6 S — 0|(p)|| ^ Po}, our choice 
of po guarantees that Uj(p) n Uj + \{p) = for p' G [oj+i — 2b,aj + \ + 26]. To define the 
extension, we introduce the following modification of the Hamiltonian N e . Let Zj{9 s ,p) 
be a function T n_1 x M n — >• R satisfying the following properties: 

• There exists C depending only on A, ||Z||c 2 and n such that ^ C. 

• Zj{9 s ,p) = Z(9 s ,p) whenever \\9 S - 9 s j {p)\\ < p . 

• Zj(9 s ,p) < Z(0 s ,p) for all 6» s and p. 

• Forpf G [aj--26,a i+ i + 26], we have that Zj(9?,pt,pf)- Zj(9 s ,pl,pf) ^ §||0 s -0j|| 2 
hold for all 9 s G TP 1-1 . 

To see that such a modification exists, let p > po be such that Zj{9 S j,p%,p^)—Zj{9 s ,p%,p^) ^ 
1110 s — 0j|| 2 on — 9 S \\ ^ p. How large p — po is depends only on A and ||Z|| C 2. Let 
Q : T 1 - 1 xl" — ► R be a smooth function such that Z(9 s ,pi,p f ) - Q = ^\\9 S - 9 s j \\ 2 for 
\\9 S - 0j|| < p and Z{9 s j ,p%,p f ) - Q ^ |||6» s - 6>||| 2 for all p/ and 6> s . The norm of Q only 
depends on A and n. Let f POj p '■ T n x R n — > R be a smooth function such that <Pp ,p = 1 
on {9 s , \\9 S — 9j\\ < po} and cp P0iP - = on {9 s , \\9 S — 9 S \\ > p}. The norm of cp POiP - depends 
only on n, po and p. Then we can choose 

Zj = (1 - ippcp) 2 " + fpo,pQ- 

We write iV £J - = H + eZj + e.R. 

For each c = p*(c^) with c/ G [aj — 26, aj+i + 26] we define 

aj (c) = aw eij (c) , .Ajv.j (c) = A Nt j (c) . 

It is not clear that these definitions are independent of the choice of the modification Zj 
or the decomposition N e = Hq + eZ + eR. We resolve these questions, and provide some 
more properties of these definition in the following proposition. 

Proposition 6.5. Let N t = Hq + eZ + eR be a Hamiltonian satisfying the genericity 
conditions [G0]-[G2] and that \\R\\ C 2 d~. There exists eo,5o > such that for < e < eo 
and < 5 < 5o the following hold. 

1. The definitions of aj and A.jsr e j(c) are independent of the decomposition N t = Hq + 
eZ + eR as long as Z satisfies [G0]-[G2] and \\R\\ C 2 sj 5; the definitions are also 
independent of the modification Zj, as long as it satisfies the same 4 bullet point 
properties. 
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2. For each c = p*(cf) with cf G [aj — 2b, aj+i + 2b], we have the local Aubry set An t ,j{c) 
is contained in the set {\\9 S — 0j|| ^ p^\ where pi is as in Theorem \5.1[ It follows 
that AN t ,j{c) C Xj and irgf \AN e ,j{c) i> s one-to-one with Lipshitz inverse. 

3. For c = p*{cf), a(c) = aj(c) if c* G [aj + b,aj+i — b]; a(c) = maxja^c), a J+ i(c)} 
if c* G (flj+i — b, aj + i + b). In particular, ay(c) > <x, + i(c) for c* = aj + ± — b and 
«j+i(c) > otj(c) for cf = Oj + i + b. 

4- For any C* G [aj+i — b, ctj + i + b], if a(c) = ay and a(c) ^ a J+ i(c), then An c { c ) = 
AN e ,j(c). Similar statement hold with j and j + 1 exchanged. 

Proof. We will prove the second statement first. The modified Hamiltonian N e> j is such 
that the single peak case of Theorem 15.11 applies, with b replaced by b/2. By choosing a 
smaller 5 if necessary, we can guarantee that p2 can be chosen the same as in Theorem 15. 11 
Theorem 15.21 also applies, where we obtain the projection property. 

We will now show that the set AN t ,j(c) depends only on the value of N € on the set 
{(8,p), \\9 S — 9 s (p)\\ ^ po}, which will imply that the definition of AN e ,j(c) is independent 
of decomposition or choice of the modification, since for all different decompositions and 
modifications, the Hamiltonian agree on this neighborhood. As before, we denote by 
N*{G,v,t) the Lagrangian corresponding to N e and JV* ■ the Lagrangian corresponding 
to N*y The projected Aubry set A]sf et j(c) is defined by the set of 9 G T n such that 
hN* .,c(G,S) = 0, where the subscript is added to stress the Lagrangian and cohomology 
class in the definition. It follows from the second statement of the proposition that any 9 
such that h]y* . )C (0 ; 9) = must be contained in {\\9 S — 9 s (c) || ^ p<i\. The following lemma 
implies independence of the local Aubry set on the decomposition or the choice of the 
modification. 

Lemma 6.1. Let N e j = Hq + eZj + eR and N e j = Hq + eZj + eR be such that 

• N e j=N e jfor \\9 S -9 s (p)\\^ Po . 

• Forpf G [a j -2b,a j+ i + 2b], we have that Zj(0 s p pl,pf)- Zj(0 s ,pi,p f ) > §||6> s -0||| 2 
and that Z^^pf) - Zj(e s ,pl,pf) > f ||0 S - 0j|| 2 for all 9 s G T n_1 . 

• ||.R||c 2 j Il-^llc* 2 ^ 

Then for sufficiently small e, 5, and for c = p*(cf) with cf G [aj — 2b, Oj+i + 2b] 

h KpC {9,9) = o h Nl . iC {e,e) = 0. 

Proof of Lemma \6.1[ Let 9q G .Ajv ?(c), we refer to Lemma [5.61 b efore and note that there 
exists an increasing sequence of integers nj., 9k = (9f,,9l) : R — > M a sequence of 
nk— periodic absolutely continuous curves such that 0^(0) = 9q and 

lim / N: i (t,9 k ,9 k )-c-9 k + a j (c)dt = 0. 

k — >ooj 
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Moreover, the curves 9k can be chosen to be minimizing, i.e. they minimizes the integral 
in the above displayed formula among the n^— periodic absolutely continuous curves such 
that 0fc(O) = Oq. In particular, 0fc|[O, rife] must be trajectories of the Euler-Lagrange flow. 
Lemma 15.61 states that for sufficiently large k we may assume that the whole curve 9kS 
are contained in \\9 S — # s (c)|| ^ po (choose a smaller 5 if necessary). 

Let (9k, Pk) be the corresponding Hamiltonian trajectory to (9k, 9k), we will show that 
for rik > \l yfi, — c\\ ^ Cy/e, where C is a constant depending only on A and n. Let 

r G [0, rifc] be where ||pfc(£)|| takes its maximum. Consider a shift 9' k (t) = 9k(t + r — l/y/e) 
of 9k, and let p' k be the corresponding action variable, then \\p' k \\ reaches maximum at 
t = l/y/e. We will write T = l/y/e in the rest of the proof. Similar to the proof of 
Proposition 15.41 we lift 9' k to a curve in W 1 without changing its name, and define 

0'k, x (t) = 9 k (x)+tx/r. 

We have the following 

[ T Nlj(t, 9' Kx , 9\ ) -c-9 k + aj (c)dt - f T N^t, 9' k ,9\) -c-9' k + a 3 (c)dt 
Jo Jo 

^ (-c + d v N: d (r, 9 k (r), 59 k (T))) ■ x + 3A^~e\x\ 2 = (-c + p' k (T)) ■ x + 3Ayfe\x\ 2 , 

the computation is identical to ([5]) and the two formulas that follows it. Assume — 
c|| > (otherwise there is nothing to prove), and we choose x to be a unit integer vector 
that minimizes x ■ (—c + p' k (T)) among unit integer vectors. We have that there exists 
C > depending on n that (-c + p' k (T)) ■ x < -C'\\p' k (T) - c||. Since 9' k (T) and 9 k (T) 
projects to the same point on the torus, by minimality of 9k we have that 

< (-c + P k(T)) -x + 3AV~e\x\ 2 < -C \\p' k (T) - c\\ + ZAsfe, 

it follows that ||p' fc (T) — c|| ^ 3A/C'y/e. Choose C = 3A/C and we have proved our claim. 

To summarize, we have proved that for sufficiently large, the curves (9k, Pk) satisfy 
ll^fc - $j( c )ll ^ P2 and \\pk~ c|| ^ C\fe. By choosing a sufficiently small e, we can guarantee 
that — 9j(pk)\\ < po- This implies that the Hamiltonians N e k and N e> k take the same 
values on the curves (9k, Pk), by taking the Legendre transform, we can conclude that the 
Lagrangian N* k and N* k must take the same values as well. It follows that 

= h Nt 3:C (9 , 9 ) = liminf / N*At, 9 k , 9' k )-c-9 k + a 3 (c)dt ^ h Re c (9 , 9 ) > 0. 

Hence h]\T e j, c (9o, Go) = =>• jC (9o, #o) = 0. The other direction also holds since 
the argument is completely symmetric. This concludes the proof of the lemma. □ 

The alpha function of a Lagrangian L can be defined by a(c) = — inf^(L — c • 9)d/i, 
where \x is taken over all invariant probability measures supported on the Aubry set A(c). 
Consider two Hamiltonians N e j and N e j as before, since the Aubry sets are identical 
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for these Hamiltonians, and the Hamiltonians coincide on a neighborhood of the Aubry 
sets, the alpha function cej(c) defined for these Hamiltonians must also be the same. This 
conclude the proof of the first statement of our proposition. 

We now prove statements 3 and 4. Consider the cohomology classes c = p*(c^) with 
c? G [a,j + b, Oj+i — b] , we note that the function Z already satisfies the conditions that 
we require of the modification, and since the local Aubry set is independent of specific 
modifications, we conclude that AnX c ) = A/v e ,j(c) and aj(c) = a(c). 

We now focus on the cohomology classes c = p*(e') with c* G [oj+i — b, Oj+i + b}. 
Using Theorem 15.11 for these cohomology classes the Aubry set An c ( c ) 1S contained in 
the vertical neighborhood {\\p — c|| ^ 36Ay/e}, and horizontally in the neighborhood 
{\\6 S -e^c)]] ^ p 2 }U{\\e s -9 s J+1 (c)\\ ^p 2 }. Takeapoint^oG^ e (c)n{||^-^(c)|| < p 2 }, 
by going through the same argument as in the proof of Lemma 16.11 we can conclude that 
^7V e ,c(#0) #o) = implies that h,N e j,c(^0j &o) = 0. It follows that AnX ) Rkj(c) C _4/v e j(c); 
the same holds for j + 1. We have that 

a(c) = — min{inf(iV* — c • 9)dfj,i,'wf(N* — c • 9)dpL 2 } ^ max{aj(c), (Xj + i(c)} 

Ml M2 

where p\ is supported on An e (c) n Uj (c) while p 2 is supported on An c (c) n Uj + i(c) . On the 
other hand, since a(c) = — m£^(N* — c ■ 6)dp with \x taken over all probability invariant 
measures, a(c) ^ aj(c), ctj + i(c). We conclude that a(c) = max{aj(c), aj + \(c)}. We have 
proved statement 3. 

Moreover, assume that An € (c) H C/j-(c) / 0, then there exists #o m this set such that 
hN etC (9o,9o) = 0, as well as a sequence of localized periodic curves 9k converging to it. 
By taking any weak-*-limit of probability measures supported on these curves, we obtain 
at least one measure v supported on An c {c) R Uj(c) such that a(c) = — J (N* — c ■ 0)du. 
This implies that a(c) ^ otj(c), hence a(c) = ctj(c). As a conclusion, if a(c) ^ otj(c) then 
An c {c)r\Uj (c) = 0. This proves statement 4 and concludes the proof of Proposition l6.5l □ 

6.2 Generic property of *Aat £ (c) 

In this section we discuss the property of the sets A/v £ (c) for c = (pKc*)) with & G 
[a,j — 2b, oti+i + 26] and their properties when we allowed to subject the Hamiltonian to 
an additional perturbation. It is convenient for us to fix a modified Hamiltonian N e j and 
base all discussions on this system. 

From Proposition 16.51 we have that the sets A/v e j(c) (we will write Aj(c) for short 
in this section) are contained in the NHIC Xj, and ngf\Aj(c) is one-to-one. We will 
study finer structures of the Aubry sets, by relating to the Aubry-Mather theory of two 
dimensional area preserving twist maps. We will prove the following statement. 

Proposition 6.6. There exists eo,<5o > such that for < e < eo and < 5 < So, there 
exists arbitrarily small C r perturbation N' e of N e , such that for each & G [aj — 2b,aj + \+2b], 
An' c {p*{c^)) supports a unique c— minimal measure. 
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We note that the time-one-map of the Hamiltonian flow is a twist map defined on 
T n x M n . The generating function of this twist map is Gj(x,x') : M. n x W 1 — > R , 

Gj(x,x')= inf f N*Jt,j,j)dt. 

7(0)=x, 7 (l)=x' J 

Consider an orbit {(9(t),p(t))} of the Hamiltonian flow, its trajectory in the configuration 
space can be lifted to a curve x(t) £ M. n , which is unique modulo integer translation. 
The sequence x^ = x(k), k £ % will be called a configuration. A configuration's rotation 

number is defined by lim^ >oo( x a+k ~ x a)/k, if such a limit exists. 

Let {x^} = {{x s k ,xj)} be a configuration corresponding to an orbit in Aj(c), we will 
say that this configuration belong to the Aubry set for short. Since Aj(c) C Xj, we have 
that the slow component x s stays bounded all the time. Take two configurations {x^} and 
{yk}, we say that they intersect in the fast direction (in short, intersect, as this is the only 
type of intersection we will consider) if there exists an integer m and indices fei, k 2 such 
that x k > yj£ + m and xj 2 < y^ + m. We have the following statements, analogous to 
the twist map case. 

Lemma 6.2. 1. Any two distinct configurations {xt} and {yk} in Aj(c) does not in- 
tersect. 

2. Any configuration {x^} in Aj(c) has a uniquely defined rotation number p = (0,pf). 

Proof. For the first statement, we prove by contradiction. Assume that x(t) and y(t) are 
the lifts of two distinct trajectories such that {x(k)} and {y(k)} intersect. It follows that 
there exists m and k\. &2 such that x* {k\) > y^(k±) + m and x^(k2) < 1/^(^2) + m. It 
follows that there exists r S M such that x* (t) = y' (r) + m. Let 9\{t) and (>2{t) be 
the projections of x(t) and y(t) to T n , we have that 0((t) = 01(t). Assume that Pi(t), 
i = 1,2 are the corresponding action variables for trajectories 9i. Let k ^ r < k + 1, we 
have that (&i(k) , Pi(k)) S Aj(c). From the graph theorem, we have that (9i(k) , Pi(k)) is a 
function of 9\ (k). Applying the flow, we have that {Oi(t),pi(t)) is a function of (9j (t),t). 
It follows that (0i(r),pi(r)) = {9 2 (t), P2 (t)), hence (0i(t),pi(t)) = (6 2 {t),p 2 (t)) for all t, 
a contradiction. 

For the second statement, since any trajectory from A/(c) must contained in Xj, 
we have that any lift x(t) of such a trajectory must have its slow component uniformly 
bounded. Hence lim^ >OQ x s (k)/k = 0. It suffices to consider only {xj,}. Since x? is one- 
dimensional, most argument from the standard Aubry-Mather theory applies, once we 
establish the non- intersecting property. We refer to \ 1 1 'j . section 11, where existence of 
rotation number was proved under a weaker assumption (the Aubry crossing lemma) . □ 

Let p be a c— minimal measure for N e j, we know that it is necessarily supported on 
Aj{c). The rotation number of p is p(p) £ i^i(T n ,M) ~ W 1 , defined by 

(c,v)dp{9,v) = (c,p(p)). 

T"xR™ 
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Using the no-intersection property (Lemma l6.2t 1), most of the statements we will be need 
follows from standard Aubry-Mather theory. Most of the arguments presented here are 
variations of those found in sec [MF] . 

Proposition 6.7. For any c = p*{cf ), & £ [aj — 2b, aj+i + 2b], the following hold. 

1. All c— minimal measures supported on Aj(c) have a common rotation number p(c) = 
(0, pf(c)). Moreover, the function a jP*(cf) as a function of c* is C 1 . 

2. If pf (c) = p/q £ Q, written in lowest terms, then all minimal measures are supported 
on q— periodic orbits. These orbits corresponds to (p,q)— periodic configurations {x k } 
in the sense that (x s k+q , %p. +q ) = i x k i x k ) + (®iP)- Furthermore, they are the minima 
of the functional 

9-1 

k=0 

over the set of configurations that are (p,q)— periodic. 

3. If p f (c) i <Q>, then there is one unique c— minimal measure. 

Proof. First we show that all the configurations on Aj(c) has the same rotation number. 
To see this, consider any two configurations with different rotation numbers, since they 
must intersect, Lemma 16.21 implies that they cannot both be contained in Aj(c). 

We now look at the function a>jp*(cf). It is known that (see, e.g. [Mai]) <x/(c) is a 
convex function and any rotation number p of a c— minimal measure is a subderivative of 
ctj at c. If for some c the subderivative is unique, then a is differentiable at c. It follows 
a j(P*( c ^)i i s differentiable for each G [aj — 2b, a J+ i + 2b]. The fact that it is C 1 
follows from the following statement: let f(x) be convex, x n is a sequence that converges 
to x*, p n is a subderivative of f{x) at x n and p n converges to p*, then is a subderivative 
of f(x) at x*. This concludes the proof of the first statement. 

We now prove the second statement. Consider any configuration {x k } with rotation 
number p/q, we have that %l +q —xl—p does not change sign for this configuration. Assume 
that it does, say Xf, i+q — x^ —p > and Xk 2+q — x k2 —p < 0, then the configurations {x k } and 
Xk+k 2 -k 1 intersects, contradiction. On the other hand, since the rotation number is p/q, 

we have that lirm. 5.00 xl +q — xl —p = 0. It follows that any x k such that xl +q — xl —p 7^ 

does not project to a point on the support of an invariant measure, since this point is not 
recurrent. By the same argument, we can show that x s k+ — x s k = for any point that 
projects to the support of an invariant measure. 

We have proved that any point on the support of an invariant measure lifts to a 
configuration with x k+q — x k = (0,p). Let p be a c— minimal measure supported on 
(9(k),p(k)), k = 0, ■ ■ ■ q — 1, and let x k be the corresponding configuration. Since 

/ (N*j - c ■ 0)dp = / N*jdp + c • P = ^2 Gj(x k ,x k+ i) + c ■ p, 
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\i minimizes J (N*j — c • 9)d/j, implies that {xk} minimizes Ylk=o Gj{xk-, %k+i)- 

For the irrational rotation number case, we refer to [MFJ, section 12. Consider Aj(c) 
as a subset of T and the dynamics on this subset. It is proved that the system is semi- 
conjugate to a rigid rotation of irrational rotation number, and the semi-conjugacy is not 
one-to-one on at most countably many points. It follows that the dynamics on Aj(c) has 
one unique invariant measure, since irrational rotation is uniquely ergodic. □ 

For irrational rotation numbers, we have that the corresponding minimal measure is 
unique. For rational rotation numbers, it is well known that for the twist map, generically 
there exists only one minimal periodic orbit of rotation number p/q. We have the same 
conclusions here. The following statement and Lemma 16.21 imply Proposition 16.61 

Proposition 6.8. 1. By subjecting the generating function Gj(x,x') to an arbitrarily 
small C r perturbation, we have that for any rational rotation number p/q, there are 
exactly q periodic configurations of type (p,q). (In this case there exists a unique 
minimal periodic orbit with rotation number p/q.) 

2. The perturbation to Gj in part 1 can be realized by an arbitrarily small C r pertur- 
bation to the Hamiltonian N e j, localized in the set {(9,p) : \\9 S — 9j(p)\\ < po}. 
As a result, this perturbation can be realized by a small perturbation to the original 
Hamiltonian N e . 

Proof. Let {xk} be a minimizing configuration of type (p,q), let U be an open set that 
contains xq but none of the x\, ■ ■ ■ , x q -\. Let gxjix) be nonnegative periodic function that 
is supported on U, gu( x o) = is the unique minimum and d 2 g is positive definite. If we 
consider the new generating function 

Gj{x,x) +gu(x), 

the action ^fc=o Gj( x ki x k+i) 1S unaffected, while the action increases for other configu- 
rations. It follows that {xfc} and its translations are the unique minimal configurations. 
However, this perturbation cannot be realized by a localized perturbation to the Hamilto- 
nian (to be more precise, it is localized horizontally, but not vertically). We consider the 
following modification of the above construction. 

Let $ denote a lift of the time-one-map of the Hamiltonian flow. The generating 
function uniquely determines the map $ in the sense that given x, x' £ W 1 , write p = 
—d\Gj and p2 = diGj then <£(x,p) = (x',p r ). On the other hand, Theorem 15.11 implies 
that any orbit in the Aubry set Aj(c) is localized in the set {||p— c|| ^ 6nA^/e}, which leads 
us to the following definition. Let V x {QA^fne) = {x' G M. n , \\d\Gj(x, x') — c\\ ^ 6^4-y/ne}, 
and let p x be a smooth function that takes value 1 on V x {QA^/ne) and takes value on 
V x (12 Ay/ne). We have that the generating function 

Gj(x,x) +gu{x)p x {x) 
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will make {xk} and its translation the unique minimizing configurations of type (p, q). The 
norm of the perturbation can be arbitrarily small since the norm of gjj can be arbitrarily 
small. 

To treat all rational rotation numbers, we consider a sequence of such perturbations 
gUi (x)p x (x'), each subsequent perturbation can be chosen to be small enough, such that the 
result of earlier perturbations are not destroyed. The final perturbed generating function 
is 

Gj(x,x') = Gj(x,x') + ^g Ui (x)p x [x'). 

We now show that the perturbation can be realized by a localized perturbation of the 
Hamiltonian. Write g{x,x') = Yli>i 9u t { x )Px{x') and let $' denote the perturbed time- 
one-map of the Hamiltonian flow. Since dig = for all x ^ \\ i Ui or x' £ V x (12A^/ne), the 
perturbed time-one-map <£' is identical to the original $ for any (x,p) ^ (Jj U x {\\p — c|| ^ 
12A^/ne}. Since we can choose Ui such that |J i Ui C {\\9 S — #j(c)|| < d < pq}, for 
sufficiently small e we can guarantee 

|J Ui x {\\p - c\\ ^ UA^Te] C {\\6 a - < p }. 

i 

It follows that $ = <&' for any (9,p) {\\6 S - 9 s j (p)\\ < p }. This perturbation of the 
time-one-map can be realized by a perturbation to the Hamiltonian localized in the same 
neighborhood. □ 

6.3 Generic property of a(c) and proof of Theorem 16.11 

After obtaining the desired properties for the local Aubry set, we now return to the 
Hamiltonian N e . If c* G [aj + 6, flj+i — b], we have that AnS c ) = •^■N € ,j{ c )- For cf G 
[oj+l — b,aj + i + b], Proposition 16.51 statement 3 and 4 shows that it suffices to identify 
whether a(c) is equal to aj(c) or Oj+i(c). 

Proposition 6.9. Assume that N t = Hq + eZ + eR is such that Z satisfy [G0]-[G2] and 
that \\R\\ C 2 ^ S. Then there exists €q,Sq > such that for < e < eo and < S < So, 
there exists an arbitrarily small perturbation N' e of N e , with the following properties. For 
the Hamiltonian N' e Proposition 1 6. 51 and Proposition \6. 6\ still hold, in addition, there exists 
only finitely many C* G [cij+i — b, oy+i + b] such that aij(p*(&), c? ) = a J+ ip*(c^). 

Proof. By taking a small perturbation if necessary, let us assume that we start with a 
Hamiltonian N e such that Proposition 16.51 and Proposition 16.61 already hold. Consider the 
interval of G [aj - 2b, a j+1 + 2b] first. Let PV{9, p, t):T n x R 71-1 x [ aj - 2b, a j+1 + 26] — ► R 
be a family of smooth functions such that 



pne, P ,t) 



'n, \\9 S - 9 S ^)\\ < po and pf G [ay - 3r/2, a i+ i + 3r/2] 
0, \\0» - 01(pf)\\ >4p /3 
k 0, p J G {ay - 26,a i+ i +26} 
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Clearly ||-Pj 7 ||c r can be arbitrarily close to by choosing 77 close to 0. 

Let N' n = N e + Pj. The new perturbation can be considered part of R and if 77 is 
sufficiently close to 0, Proposition 16.51 still hold. This implies that the local Aubry sets 
still depends only on the value of the Hamiltonian on the set \\9 S — 0j(pf ) || ^ po, on which 
the perturbation is simply a constant (for & G [dj — 3r/2, Oj+i + 3r/2]). We have that 
aNvj(c) = ctNt,j + V and that Anvj(c) = -4jv e y(c) for c* G [oy — 3r/2,ay + i + 3r/2]. It 
follows that all properties of the local Aubry set AN e ,j(c) is intact, while ctj(c) undergoes 
a shift. 

On the other hand, Consider the functions ajp*(cf) and aj+ip*(c-f) as functions on 
[cLj + \ — 3r/2, Oj + i + 3r/2]. Since they are both C 1 , by Sard's lemma, the critical values of 
ay — ay+i has zero measure. It follows that there exists a full measure set of 77 G K such 
that a'j — a'- +1 = implies otj — ay+i + rj ^ 0. In other words, the two functions ay + and 
ay + i intersect transversally, which implies that there are only finitely many values where 
ay - Oj+i + 77 = 0. 

We can perform this modification for each [a,- — 2b, oy+i + 26], and 77 can be chosen to 
be arbitrarily close to 0. □ 

Proof of Theorem \6.1[ Since there are only finitely many G [dj+i — b,aj + \ + b] on 
which otj = aj+i, we add these points to the set {ao,--- ,a s } to form a new partition 
{[<%, On each open interval (cij, Oj+i) a(c) is only equal to one of the ctj and a J+ i. 

Use Proposition 16.51 and the first statement follows. 

The second statement follows from Proposition 16.61 □ 

6.4 nondegeneracy of the barrier functions 

In this section we prove Theorem 16.41 We have concluded that in order to prove Theo- 
rem [631 it suffices to show that Ti(e) = r*(e) and Tz^e) = T^e). We show that this is 
the case by proving the following equivalent statement. 

Proposition 6.10. Let H' e be a perturbation of H t such that the conclusions of Theo- 
rem \6.2\ holds, then there exists an arbitrarily small C r perturbation H" to H' e , such that 
for the Hamiltonian H" Theorem 1 6. 3 still hold, in addition: 

1. Consider c$ £ (cij, such that A c = N c and TtgfAc = T. Take ( G Ai c , and let Ci 
and (2 be its lifts to the double cover. We have that the functions h c ((i, 9) + h c (6, C2) 
and h c (C,2,6) + h c (9,C,i) have isolated minima outside of the lifts of A c . 

2. For c = Cij+i, take ( G A c n Xj and rj G A c H Xj + \. We have that both h c (C,9) + 
h c (6,rj) and h c (rj,9) + h c (6,() has isolated minima outside of A c . 

This proposition is essentially proved by Cheng and Yan in |CY2j . here we briefly 
describe their approach. 

Consider the Hamiltonian H' e such that conclusions of Theorem 16.21 holds. In the 
rest of the section, let's refer to H' e simply as H. For now, let us also fix an interval 
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(a,j, fflj+i) and consider only cohomology classes with cf in that interval. Let = Ti(e) n 
{c,c-f 6 (cbj, aj+i)}, we would like to show by perturbing the Hamiltonian, we can make 
the functions h c (£i,9) + h c (9,^2) and h c (^2,9) + h c (0,£i) nondegenerate. 

Recall that h c is the barrier function defined on the covering space (2T) n x R n , and 
£ : (2T) n x M. n — > T n x R™ is the covering map. H is the Hamiltonian lifted to the 
covering space. 

Define the generating function G(x, x') : W 1 x R n — > R by 

G(x,x')= inf / 1 L(t, 7 ,7), 

7(0)=x,7(l)=x' Jo 

where L is the Lagrangian corresponding to H. A convenient way of introducing per- 
turbations to the functions h c is by perturbing the generating functions. Denote by 
7r : R™ — > T n the standard projection. 
We consider the following perturbation 

G'{x,x') = G(x,x , ) + Gi(x / ) 

and denote by h' c the corresponding perturbed barrier function. We have the following 
statement. 

Lemma 6.3. flCYZtf . Lemma 7.1) For c = p*(c^) with cf G (a,j, %+i), the following hold. 

1. There exists a family of open sets O c C (2T) n such that the full orbit of any (9,p) £ 
A/&(c) \ »4&(c) musi intersect O c in the 9 component. 

2. There exists p > such that if we perturb G(x, x') by function G\ (a/) with supp G\ C 
B p (u), where B p (u) is the ball of radius p centered at u, then for each c such that 
Bp[u) C 7r -1 O c the corresponding barrier function 

h'Mi.O) + ~h' c (9, 6) = hcituO) + hc(0, 6) + <n (0) 

/or eac/i 9 G O c . 

5. £O c n {6» : ||<9 S - 0|(c)|| < p 2 } = 0, in particular, £O c nJV fl (c) = 0. Moreover 
# = U c / e ( aj ,a J+1 ) °c ^ an open set. 

As before, let us write b+(9) = h e (£ 1} 6) + h c (9,&) and b~(9) = h e (&,0) + h c {9,^). 
Elements of A/fv(c) \ «4fr(c) coincide with the minimal set of the functions b^r. To prove 
that this set is isolated, it suffices to prove its intersection with O c is isolated, as any 
accumulation point of Njj[c) \ Ajj(c) has a diffeomorphic image in O c . We say that the 
function b^(9) is degenerate if its minimal set has at least one accumulation point. Cheng 
and Yan proved that it is possible to introduce a perturbation to make b^r nondegenerate 
for all c G rj simultaneously. 
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This is not possible in general, if the functions b^r behave badly as c varies. Since 
regularity of b^r in c is hard to prove, Cheng and Yan resolves this problem nicely by 
introducing an additional parameter. Recall that for each c G T\, the Aubry set A is an 
invariant curve on the time-zero section of the invariant cylinder Xj, call it 7 C . Fix an 
arbitrary curve 70 = {pf = Pq} H {t = 0} n Xj, we introduce a parameter a which is the 
area between 7 C and 70 on the cylinder Xj n {t = 0}. a is monotone in (/ and is only 
defined for c £ Tj. Cheng and Yan proved that 

Lemma 6.4. ([CY2], Lemma 6.4) There exists constant C > such that, for a and a' such 
thatc(a),c(a')er{,CeA H (c) and m £ {\\0 S - 0<(c(ff))|| < p 2 }U{||0 s -0j(c(<7'))|| < p 2 }, 

\h c (a){(,m) ~ K(a')(C,m)\ «S C(\/k - + l c (°") - c^OD) 

|^c(o-)( m )C) - ^ c (<T')( m >OI ^ C(y/\cr - a'\ + |c(tr) - c(cr')l). 

It follows that the function h c ^ can be extended to /i Cj(J that is ^—Holder in c and cr, 
this regularity turns out to be enough. To see how this is carried out, let us consider a 
subset Bd'(c[) x Rd(u), where Bd>{cl) = {c : \cf — c{\ < d'} and Rd(u) C T n is an open 
cube centered at u with edge d. 

Lemma 6.5. ([CY2], Lemma 7.2) There is a residue set of functions G\ G Cq(R<i(u), H.) 
such that 

bt(e) + G 1 (e) 

has isolated minima in Rd(u) for each c G Y\ n Bd'(c). (Cq stands for C r functions with 
compact support). 

Remark 6.1. The nontrivial part of this statement is that the nondegeneracy of can 
be achieved for all c G T\ n B^ (c) simultaneously. The regularity acquired in Lemma \6.4\ 
is crucial to the proof. We refer to JCY2] for details. 

To construct the desired perturbation to the barrier function, let us state another 
lemma, which is a consequence of the upper semi-continuity of the Mane set on the La- 
grangian. 

Lemma 6.6. The property that the functions b^r are non- degenerate on the set Bd'(cl) x 
Rd(u) survives under sufficiently small perturbation. 

We proceed to prove Proposition IfTTUl Let B d ,{c{) x R di (ui) C U, be a sequence 

of sets such that U = {JiB d '_(c{) x We may choose a sequence of perturbations 

Gi : Rdi (u{) — > R, and let G' k (x, x') = G(x, a:')+X^=i Gi( x ') an d bf k be the corresponding 
barrier functions corresponding to the generating function G' k . We can choose the sequence 

Gi inductively such that bf is non-degenerate on (c, 9) G Ui=i(-^d' ( c { ) ^ ^1) x Rdi( u i), 
because new perturbations can be added that does not disrupt the nondegeneracy already 
established in the previous steps. By repeat this process for each interval (a,j,dj + i), we 



48 



have constructed a perturbation to the generating function G, such that the first statement 
of Proposition 16.101 holds. 

For the second statement, using the same arguments for Lemma 16.31 one can show 
that the same type of conclusions apply to b^r as well. 

Lemma 6.7. For each c = p*(c*) with = cLj,j = 2,--- ,1 — 1 the following hold. 

1. There exists a family of open sets O c C (2T) n such that the full orbit of any (9,p) G 
Mh{c) \ Ah{c) must intersect O c in the 9 component. 

2. There exists b > such that if we perturb G(x,x') by function G\{x') with suppGi C 
Bb(u), where Bb(u) is the ball of radius b centered at u, then for each c such that 
Bfy(u) C 7r _1 O c the corresponding barrier function 

ti c ((, 9) + h' c (9, rj) = h c (C, 9) + h c (9, n) + G 1 (9) 

for each 9 G O c . The same conclusion holds for h c (r],9) + h c {9,C). 

For a fixed c, it is easy to see bf(9) + G\{9) has isolated minimal set in Rd{u) for an 
open and dense set of G\. Repeat the arguments for bf , we obtain a perturbation for 
which the both statements of Proposition I6.1UI hold. □ 

A Generic conditions 

We prove Theorem 12. II in this section. Consider the following (degeneracy) conditions on 
the function Z(9 s ,p*(pf)) : T n_1 x [a min , a max ] — > R. 

[TO] For pf £ [a 

mini ftmax] > &11 local maxima of Z {9 s , p*(pJ ')) is nondegenerate. 

[Tl] For each pf £ [a m i n , a max ] and there are at most two distinct 9{ , #| £ T n_1 such that 
d e sZ{9 s j ,p^{pf)) = for j = 1,2 and that Z{9 s 1 ,p^{pf)) = Z{9 s 2 ,p*{pf)). 

[T2] For any pf G [a min ,a max ] and distinct 0?,0| G T n_1 such that d 8 sZ(9 s j ,p^(pf)) = 
for j = 1,2, we have that 

d p fZ{9lp*{pf)) ^ d pf Z{9 s 2 ,p*{pf)). 

Let W denote the set of functions in S r that satisfies one or more of the conditions [TO]- 
[T2]. 

Proposition A.l. W is open and dense. 

Proof of Theorem \2.1\ The set IA is open, since if H\ satisfies conditions [G0]-[G2] with 
some A > 0, any H[ sufficiently close to Hi in C r norm satisfies these conditions with a 
slightly smaller A' > 0. 
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We now prove that U is dense by showing that U DU'. The conditions [T0]-[T2] implies 
the statement that any pf is either a nondegenerate regular point or a nondegenerate 
bifurcation point, and that there are at most finitely many bifurcation points. To see that 
[T0]-[T2] also imply conditions [G0]-[G2], let {[dj, Q>j+i]}jZo be the partition of [a m i n , a max ] 
by bifurcation points. Each pf G (aj,a,j+i) defines a unique global maximum 9 s Ap f ). The 

function 6 s {pf) is continuous since any converging sequence 9 s (pf.) also converges to a 
global maximum, and it must be smooth by implicit function theorem. The function 
extends to [a^Oj+i] by continuity, and using the nondegeneracy of the maximum and 
implicit function theorem, we can extend 9? smoothly to the interval [ctj — d,aj + \ + d], 
such that each 9j(pf) is a nondegenerate local maxima. Assume that for each pf G 
[cij — d, aj + i + d] we have —dg sgs Z(9 s ,p*(pf)) ^ d'l as a quadratic form, hence Z satisfies 
[GO] with A = mm{d,d'}. [Gl] and [G2] are direct consequences of [T0]-[T2]. □ 



B Normally hyperbolic manifold 



Let F : W 1 — > R n be a C 1 vector field. We give sufficient conditions for the existence of 
a Normally hyperbolic invariant graph of F. We split the space R n as R ra " x R" s x R nc , 
and denote by x = (u, s, c) the points of R n . We denote by (F u , F s , F c ) the components of 
F: 

F(x) = (F u {x),F s {x),F c {x)). 
We study the flow of F in the domain 

Q = B U xB s x Vl c 



where B u and B s are the open Euclidean balls of radius r u and r s in W""" 
£l c is a convex open subset of M. nc . We denote by 



and M™ s , and 



L(x) = dF(x) 



! uu(x) L us {x) L uc {x) 



J-'suix) L ss (x) L sc (x 
L cu (x) L cs (x) L cc (. 



x 



the linearized vector field at point x. We assume that ||L(x)|| is bounded on f2, which 
implies that each trajectory of F is defined until it leaves ft. We denote by W c the union 
of full orbits contained in Q. In other words, this is the set of initial conditions x G fl 
such that there exists a solution x(t) : R — > Q of the equation x = F(x) satisfying 
x(0) = 0. We denote by W sc the set of points whose positive orbit remains inside fi. In 
other words, this is the set of initial conditions x G fl such that there exists a solution 
x(t) : [0, oo ) — > fi of the equation x = F(x) satisfying x(0) = 0. Finally, we denote by 
W uc the set of points whose negative orbit remains inside fl. In other words, this is the 
set of initial conditions x G ft such that there exists a solution x{t) : (oo,0] — > ft of 
the equation x = F(x) satisfying x(0) = 0. These sets have specific features under the 
following assumptions: 
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Hypothesis B.l (Isolating block). We have: 

• F c = onB u xB s x dn c . 

• F u (u, s,c)-u>0 on dB u x B s x fi c . 

• F s (u,s,c) • s < on B u x dB s x C . 

Hypothesis B.2. There exist positive constants a, m and M such that: 

• L uu (x) ^ al, L ss (x) ^ —al for each x £ £1 in the sense of quadratic forms. 

• ||L us (x)|| + ||L uc (x)|| + ||L su (x)|| + ||L sc (a;)|| + ||L cu (a;)|| + ||£ca(a:)|| + ||i cc (^)|| < m 
for each x G f2. 

Theorem B.l. Assume that Hypotheses I B. 1\ and \B.S\ hold, and that 

m 1 

K := sC — . 

a — 2m s/2 

Then the set W sc is the graph of a C 1 function 

w sc :B s xff^ B u , 

the set W uc is the graph of a C 1 function 

w uc : B u xQ c — ► B s , 

and the set W c is the graph of a C 1 function 

w c = «, w c s ) : n c — »• B u x B s . 

Moreover, we have the estimates 

\\dw sc \\ < K, \\dw uc \\ ^ K, \\dw c \\ s: 2K. 

Proof. This results could be reduced to several already existing ones, see [El IHPSt IMcG} 
ICh] or proved directly by well-known methods. We shall use Theorem 1.1 in |Ya] which 
is the closest to our needs because it is expressed in terms of vector fields. We first derive 
some conclusions from the isolating block conditions. We denote by ir sc the projection 
(u, s, c) i — > (s, c), and so on. 

Lemma B.l. If Hypothesis \BJ\ holds, then 

tt uc (W uc ) = B u x n c and ir sc (W sc ) = B s x Q c . 
Moreover, the closures of W sc and W uc satisfy 

W sc C B s x B c x fl c , W uc C B s x B° x Cl c . 
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Proof. Let us define T + (x) G [0, oo] as the first positive time where the orbit of x hits the 
boundary d£l. Let us denote by ip(t,x) the flow of F. If T(x) < oo, we have (p(T(x),x) 6 
dB u x B s x Q, as follows from Hvpothesis lB.il Then, it is easy to check that the function 
T is continuous, and even C , at x. 

We prove the first equality of the Lemma by contradiction, and assume that there 
exists a point (u, c) 6 i?" x Q c such that VF" C does not intersect the disc {u} x B s x {c}. 
Then, the first exit map 

B s 9s^ ^(T(x),x) 6 55 s 

extends by continuity to a continuous retraction from B s to its boundary dB s . Such a 
retraction does not exist. The proof of the other equality is similar. 
Finally, we have 

W uc C B u x B s x n c = (B u x B s x n c ) (J (dB u x B s x Cl c ). 

Hypothesis IB. II implies that each point of dB u x B s x Q c has a neighborhood formed of 
points which leave f2 after a small time. As a consequence, the set dB u x B s x Cl c can't 
intersect W uc , and we have proved that W uc C B u x B s x £l c . The other inclusion can be 
proved in a similar way. □ 
In order to prove the statement of the Theorem concerning W sc , we apply Theorem 
1.1 of [Yaj . More precisely, using the notation of that paper, we set 



a = u/K, z = (s,c), f(a,z) = F u (Ka,z)/K, g(a, z) = {F s (Ka, z), F c (Ka, z)). 
We have the estimates 



L.ss F sc 
Fes L cc 



< m 



d a f = F uu ^ a, d z g 

in the sense of quadratic forms. Moreover, we have the estimates 

vn 

\\d z f\\^-, \\d a g\\^Km. 

Since 

m + m/K + Km < 2m + m/K = a 

we conclude that Hypothesis 2 of [Yaj is satisfied. Hypothesis 1 of |Yaj is verified by the 
domain f2, and Hypothesis 3 is precisely the conclusion of Lemma IB. II As a consequence, 
we can apply Theorem 1.1 of [Yaj . and conclude that the set W sc is the graph of a C 1 
and 1-Lipschitz map above B s x Q c in (a, z) coordinates, and therefore the graph of a 
-fT-Lipschitz C 1 map w sc : B s x Q c — )■ B u in (u, s, c) coordinates. 

In order to prove the statement concerning W uc , we apply Theorem 1.1 of [Ya] with 

a = s/K, z = (u,c), 
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f(a,z) = -F s (Ka,z)/K, g(a,z) = -(F u (Ka,z),F c (Ka,z)). 

It is easy to check as above that all hypotheses are satisfied. 

Let us now study the set W c = W sc n W uc . First, let us prove that W c is a C 1 graph 
above Q c . We know that W sc is the graph of a i'T-Lipshitz C 1 function w sc (s,c) and that 
W uc is the graph of a i^-Lipshitz C 1 function w uc (u, c). The point (u, s, c) belongs to W° 
if and only if 

u = w sc (s, c) and s = w uc (u, c), 
or in other words if and only if (u, s) is a fixed point of the i<C-Lipschitz C 1 map 

(u,s)^(w sc (s,c),w uc (u,c)). 

For each c, this contracting map has a unique fixed point in B u x B s , which corresponds to 
a point of W sc n W uc . It follows from Lemma [B , 1 1 that this point is contained in B u x B s . 
Then, it depends in a C 1 way of the parameter c. We have proved that W c is the graph of 
a C 1 function u; c . In order to estimate the Lipschitz constant of this graph, we consider 
two points (ui, Si, Ci),i = 0, 1 in T. We have 

IN - Woll 2 < ^ 2 (ll«l - Soil 2 + || C X - C || 2 ) 

and 

\\si - soil 2 < K 2 (\\ui - u \\ 2 + ||ci - c || 2 ). 

Taking the sum gives 

(1 - K 2 )(\\ Ul - n || 2 + ||*i - soil 2 ) ^ 2iT 2 ||c 1 - c || 2 

and 

/ 2K 2 

||(ui,si) - (u ,so)|| < A/ i _ k 2 ~ C °" ^ 2K W Cl ~ C °H' 

since K ^ 1/V2. We conclude that w° is 2X-Lipschitz. □ 
We need an addendum for applications: 

Proposition B.2. Assume in addition that there exists a translation g o/M™ c such that 

g(Sl c ) = ft c and F o (id <g> id ® g) = F. 

Then we have 

w sc o (id (g> g) = w sc , w uc o (id®g) = w uc , w c o g = w c . 



Proof. It follows immediately from the definition of the sets W sc , W uc and W c that 
g(W sc ) = W sc , g(W uc ) = W uc and g(W c ) = W c . □ 
In applications the first condition of Hypothesis IB, II is usually not satisfied, except in 
the case where f2 c = W lc . It is thus useful to state a more "applicable" variant of the 
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result. In view of the applications we have in mind, it is useful to split the central variables 
into two groups and consider the case 

n c = R n ' x n C2 , 

where f2 C2 is a convex open set in R n ^, nl + v? = n c . Given a positive parameter a, let 
OS 2 be the set of points c 2 € R nc such that d(c,Q C2 ) < a. This is a convex open subset 

2 1 

of W 1 ^ containing £l c ' 2 . We denote by the product K™ c x Q, C2 and by Q, a the product 
B u x B s x VL%. With the notation F c = (F C1 ,F C2 ), we have: 

Proposition B.3. Assume that there exists A, m, a > suc/i t/ioi 

• F u (u,s,c) ■ u > on dB u x B s x n c a . 

• F s (u,s,c) ■ s < on B u x dB s x n c a . 

• L uu (x) 5? a/, L ss (x) sj —a/ /or eac/i x G f^o- ift fie sense of quadratic forms. 

• || -|-||Z.^ C (^) || H-||^ ss (ic) || -|-|| J^ sc (a7) || ||Z. c ^(sc) || H-||Z. CS (a?) 1 1 + 1 1 Z. cc (a?) 1 1 -(-2 1 1 ^ C2 (a?) || /cr ^ 
m /or eac/i x G CT . 

Assume furthermore that 

TTh 1 

if := < -=. 

a — 2m ^2 

T/ien i/iere exists a C 1 function p : Q£ — >• [0, 1] which is equal to 1 on f2 c ane? smc/i i/iai 
the vector field 

F(u, s, c) := (F u (u, s, c 1} c 2 ), F s (u, s, c x , c 2 ), F C1 (u, s, ci, c 2 ), p(c 2 )F C2 (u, s, c x , c 2 )) 
satisfies all the hypotheses of Theorem \B.l\ on iVoie ftai F = F on Q. 
Proof. We take a function /> : $7£ 2 — >• [0, 1] such that : 

• p = near the boundary of f2£ 2 , 

• p = 1 on fl c " , 

• || dp || ^ 2/cr uniformly. 

Denoting by L** the variational matrix associated to F, we see that 

L cu (u, s, c) = p(c 2 )L cu (u, s, c), L cs (n, s, c) = p(c 2 )L cs (u, s, c), 
L ClCl (u, s,c) = p(c 2 )L ClCl (u,s,c), L ClC2 (u,s,c) = p(c 2 )L ClC2 (u, s,c), 

and 

L C2C2 (u,s,c) = p(c 2 )L C2C2 (u,s,c) + dp(c 2 ) ® F C2 (u,s,c). 
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As a consequence, we have 

\\L us {x)\\ + ||Z uc (x)|| + ||L„(z)|| + ||L«.(a;)|| + \\L cu {x)\\ + \\L cs {x)\\ + ||L cc (x)|| 
=p(c 2 ){\\L us {x)\\ + \\L uc {x)\\ + \\L ss (x)\\ + ||L sc (x)||||L cu (x)|| + ||L„(x)|| + \\L cc (x)\\) 
+ \\F C2 (x)\\\\dp(c 2 )\\ ^ m. 

□ 

Under the hypotheses of Proposition IB, 31 the sets W sc , W uc , W c associated to F are 
graphs of C 1 functions 

w sc :B s xn c r ^ B u , w uc :fi"xS]^ B s , w c : Sl c T — > B u x B s 

which satisfying the estimates 

\\dw sc \\ ^ K, \\dw uc \\ sC K, \\dw c \\ s: 2K. 

The restrictions to £1 

W sc _ y^sc p yyuc _ y^uc R ^ yyc = yysc R yyuc _ yyc p ^ 

are weakly invariant by F in the sense that this vector field is tangent to them. They 
satisfy various interesting properties. For example, each F-invariant set contained in f2 is 
contained in W c . 
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